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ABSTRACT 


Retailing is becoming an increasingly important area of management attention and 
academic research. The development of decision support models to help retailers 
improve their decision-making is gaining importance. Every week motion picture 
exhibitors have to make an important decision regarding the replacement of the 
movies playing at the screens in their theaters. The dynamic and uncertain decision 
environment, complicated contract terms, a number of new products and their 
perishing demand give rise to the complexity of the problem. In several instances, this 
complexity is further increased by availability and choice of mupltiple display 
facilities. 

In this thesis, we develop modelling and heuristic methods to help the 
managers in this industry improve their decision making regarding the optimal use of 
their display space. A Markov decision process based model is used to address the 
problem of exhibitors with single screen. We propose a Genetic Algorithm based 
heuristic to address the similar problem with multiple screens. We investigate the 
improvement shown by our heuristic over other heuristics. We also discuss the 
extensions of an integer programming model SilverScreener, to address similar 
problems in other application areas. 


(iv) 



CHAPTER 1 


INTRODUCllON 

Marketing is “the process of planning and executing the conception, pricing, 
promotion, and distribution of ideas, goods and services to create exchanges that 
satisfy individual and organizational objectives.” (Kotler, 2000, P. 8). Often marketing 
systems are quite complex so that marketing researchers need to develop analytical 
models to manage these systems. 'I'here are two basic methodologies for modeling in 
marketing: verbal and mathematical. Most of the models in the so-called “behavioural 
literature” in marketing are cast as verbal models. For example, Howard and Sheth’s 
(1969) theory of consumer behaviour is a verbal model of consumer behaviour. 
Mathematical models, popularly known as analytical or quantitative models, use 
symbols to denote marketing variables and express their relationship as equations or 
inequalities. The analysis follows the rules of mathematical logic. For example, 
Bass's (1969) model of product diffusion is one of the most frequently used 
quantitative models in marketing. 

Quantitative modeling in the field of marketing is progressing rapidly. These 
promising analytic and databased findings are important to the development of a 
marketing theory that enables marketers to understand and manage markets more 
effectively in an increasingly dynamic and complex environment. In particular, there 
has been a rapid development in the research addressing the issues related to dynamic 
decision making. Examples can be found in areas such as advertising (Erickson 1995), 
brand choice (Erdem and Keane 1996), new product development (Bayus 1995), 
channel coordination (Chintagunta and Jain 1992), consumer behavior (Krishna 

1994) , market research (Bockholt and Dillon 1997), pricing (Raman and Chatterjee 

1995) , promotion (Papatla and Krishnamurthi 1996), and retailing (Krider and 
Weinberg 1997). 

Retailing, one of the important functions of marketing, is becoming an 
increasingly important area of management attention and academic research, 
especially in the developed countries. One of the major areas of research in retailing is 
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the development of decision support models to help retailers improve their decision- 
making in the dynamic and complex environmenl. Most previous marketing research 
in this regard has been concerned with consumer durable goods, not with services or 
intangibles. This thesis seeks to advance our understanding of how quantitative model 
building can improve marketing decision making in complex dynamic environments 
by focusing on a perishable entertainment product, namely, movies. 'I'lioiigh we focus 
on movies in this thesis, our methodology and results can readily be generalized to 
other entertainment products (e.g., performing arts, books, video games), travel 
.services, fashion goods, and educational programs. 

Shelf (display) space management is an area of important concern to retailers. 
First, consumers choose from the products that are displayed on the shelf. In this 
sense, shelf space management determines the ultimate profitability of the retailer. 
Second, most retailers have limited shelf space to display their products. Therefore, 
the choice of which products to display becomes an important retailing decision. The 
decision is particularly complex in those industries in which many new products are 
continually introduced in market. Finally, anticipating and adapting to dynamic 
changes in consumers’ tastes and demand are key concerns to most retailers. 
Moreover, such decisions have to be made for several product categories. Though 
previous researchers (Corstjens and Doyle 1983; Bultez and Naert 1988; Borin, 
Farris, and Freeland 1994) have addressed some of these issues, they have generally 
been in the context of packaged goods in a supermarket chain setting. However, other 
industry settings, such as the motion picture (or movie) industry considered in this 
thesis, pose different challenges and intriguing problems. 

A movie is an interesting product for several reasons. Most movies are 
separate entities and can be considered an innovation because of the new features 
(e.g., actors, storyline, music, etc.) usually included in them. These new movies are 
released every week throughout the year. The product life cycle of movies is 
relatively short and is measured in weeks. The product is seasonal in nature and the 
prominent seasons are during holidays. Finally, one of the most interesting attribute of 
movies from a research standpoint, and a major component of the complexity 
associated with the problems addressed in this thesis, is their perishability. 



In response to the dynamics and challenges posed by the above characteristics 
of movies, a stream of research, particularly addressing the marketing t>f movies has 
contributed to the marketing literature. At the consumer behavior level, some of the 
research has questioned the relevance of the traditional information-seeking 
framework for studying the consumption of movies (e.g., Hirschntan and Holbrook, 
1982; Holbrook and Hirschman 1982): Another stream of research has focused on 
forecasting the enjoyment of movies at the individual level (Eliashberg and Sawhney 
1994) as well as forecasting the commercial success of movies at the aggregate level 
(Smith and Smith 1986; Austin and (Jordon 1987; Dodds and Holbrook 1988; 
Sawhney and Eliashberg 1996; Eliashberg and Shugan 1997). Additionally, some 
research has begun to emerge addressing diflusion (Mahajan, Muller, and Kerin 1984; 
Jones and Ritz 1991), seasonality (Radas and Shugan 1995), release timing (Krider 
and Weinberg 1998), clustering (Jedidi, Krider, and Weinberg 1998), sequential 
products (Lehmann and Weinberg 1998; Prasad, Mahajan, and Bronnenberg 1998), 
contract design (Swami, Lee, and Weinberg 1998), scheduling (Swami, Eliashberg 
and Weinberg 1999; Eliashberg, Swami, Weinberg and Wierenga 2001) and the 
impact of advertising (Zufryden 1996). This thesis advances the above stream of 
research in marketing of movies, by considering the problem related to retailing of 
movies. 

Typically, perishability is thought of in terms of physical deterioration of a 
product such as a grocery item with an expiry date. In some sense, this view comes 
from the supply side of the product. In contrast, in this thesis we adopt a “demand side 
view" in the context of perishability. Thus, the physical product in the problems 
considered remains the same, but its demand perishes over time. In recent years 
considerable work has been done on the treatment of perishability in inventory control 
(e.g., Abad 1996, Jain and Silver 1994). In inventory control, perishability refers to 
the physical deterioration of units of a product. For the movies, we analogously define 
perishability as the decrease in the value/appeal of a movie with the passage of time. 
As mentioned earlier, previous research in shelf space management and dynamic 
decision making has mainly focused on non-perishable goods. In the context of 
movies, the complexity introduced by perishability in dynamic shelf space 
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management problem can be summarized as; which motion pictures to choose to 
show each week and for how long to play them. In case of exhibitors with multiple 
screens, issues such as allocation of different movies to different capacity screens, 
switching of the movies between screens and multiple screening of same movie add to 
the complexity of the general problem. This decision-making problem is further 
complicated by additiimal factors specific to the movies context. Wo discuss such 
factors in detail in the later chapters to follow. 

The exhibitors, w'ho are the retailers in movie market, also face the problem of 
limited “shelf-space” which is tlie available screens for playing movies. 'I'hc problem 
has two parameters; the number ol" screens available with the exhibitor, either one or 
multiple, and the nature of perishability of demand of the movies, which can either be 
stochastic or detenninistic. Detenninistic dynamic programming, or an integer 
program, readily models a single screen deterministic demand problem. The 
incorporation of stochasticity, however, requires more specialized techniques. 

In the proposed research we address exhibitor’s problem in two ways. We 
propose an application of “Markov decision process” (MDP) assuming a single screen 
with the exhibitor and considering stochasticity in the demand of movies. The 
complexity in this problem comes from the stochasticity in the demand of movies and 
the perishability of this product in terms of declining demand. To test the proposed 
approach we compare its performance with some other heuristics. We also investigate 
the optimal replacement policies under different simulated scenarios and characterize 
the form of policies that emerge. We explore the possibility that the optimal policy 
possesses a simple structure in movie replacement problem. 

Next w'e address the problem of exhibitor with multiple screens. We assume 
an exponentially declining demand pattern of each movie for this multiple-screen 
problem. The complexity in such problems is because of multiple unequal capacity 
screens and the decay of movies which causes the exponentially declining demand 
pattern. With an assumption of equal screen capacities, the problem is similar to the 
conventional parallel machine scheduling problems. Some work has been done in this 
regard particularly focusing on the movie exhibitors (Swami, Eliashberg, and 
Weinberg 1999). However this assumption of equal screen capacities poses serious 
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limitation in implementation of such models. Considering different capacities for 
screens introduces nonlinearity into the resulting model and makes it less tractable 
(Saxena, 2000). Wc propose the use of heuristic based on Genetic Algorithms 
methodology (Holland 1975, Michalewicz 1992, Grefenstette 1986). We also 
investigate the improvement of proposed GA-based heuristic over other heuristics. 

The rest of the thesis is organized as follows. Chapter 2 reviews the relevant 
literature. Chapter 3 explains the model for using MDP in the single-screen- 
stochastic-demand problem. We test the model for various randomly generated 
scenarios, by use of design of experiments, regression analysis and ANOVA. For 
these scenarios we also compare the performance of this model, with generally 
accepted heuristics. Finally we present the implementation of this model by using 
actual data at a theatre. 

Chapter 4 discusses the SilverScreener model proposed by Swami, Eliashberg, 
and Weinberg (1999) to address the problem of exhibitors with multiple screens with 
assumptions of deterministic demand pattern of movies and equal capacity of all 
screens. Saxena (2000) explains the implementation of this model and proposes a 
conceptual nonlinear model to address this problem with an added constraint of 
unequal screen capacities. Saxena (2000) also proposes a heuristic based on demand 
of individual movies, to address the problem of unequal screen capacities. We explain 
the proposed GA based heuristic and compare its performance against Saxena’s 
(2000) Screen Allotment heuristic. We test the model for randomly generated 
scenarios and carry out ANOVA to study the impact of various factors on the 
improvement achieved by proposed heuristics over Saxena’s (2000) heuristic. 

In Chapter 5, we discuss the conceptual extensions of the SilverScreener 
model to address some problems in similar areas and propose the use of GA-based 
heuristics for problems that cannot be addressed by simple extensions of 
SilverScreener model. We conclude in Chapter 6 by discussing limitations of the 
current research. 
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We now present the summary of two main chapters contained in this thesis. 


Chapter 3: Movie Replacement Problem; MDP Approach 

The chapter deals with a shelf space management problem, with an objective of 
maximizing retailer revenue by optimum use of his/her limited display space. The 
retailer at each decision epoch makes a decision of which product to display. The 
chapter assumes a single shelf for display with the retailer, which can be allotted to 
only one product at a time. Each product when displayed will lead to probabilistic 
demand, fhe demand for other products will also deteriorate probabilistically with 
each decision epoch. I'he stochasticity in the demand for each product induces 
complexity in the problem of space allocation. 

We model this problem as a ‘Markov Decision Process’ model (Puterman 1994). We 
define a rank of a product as the expected revenue that it will produce if displayed. 
We categorize products into groups. Each group has a probability matrix that defines 
the probability of the product opening into a particular rank and the probability of 
transition of product from one rank to another, alter each stage. Eor the sake of 
simplicity we assume these probabilities of transition as static, that is, they do not 
change with the decision epoch. We also assume independence of these transitions for 
a product with respect to other product. 'I'his means that the transition of a product 
from one rank to another rank is independent of the transition of other products. 

We focus our study on motion picture industry. We investigate optimal replacement 
policies as suggested by the model under different simulated scenarios for their 
managerial implications. Next, we characterize the form of the optimal policies that 
emerge. In particular, we explore the possibility that the optimal policy possesses a 
simple structure in the movie replacement problem. We also compare the performance 
of our model with the generally accepted heuristics. 

Chapter 4: Multiple Screen Problem: Genetic Algorithm Approach 
The analogy of movie screen allocation with the shelf space allocation problem is 
further explored in this chapter. In practice, retailers allocate space to multiple 
products rather than a single product. This is also applicable to most movie theaters in 



major cities of the world that are multiplexes, that is, theaters with multiple number of 
screens. The retailers in such situations have to make two decisions. First, which 
products should be displayed, and second, how much space to be allotted to each 
product. The complexity in this problem comes from the multiple display spaces 
available with the retailer. I'he movie exhibitors have multiple screens with unequal 
capacities. For the sake of simplicity, we assume a deterministic demand pattern (with 
random error term) of movies in this case. I'hc number of visitors for a movie 
scheduled at a particular screen in a week is limited to the minimum of the screen 
capacity and the general demand for the movie in that week. This situation is 
analogous to any other retail space allocation problem. With a particular revenue 
function, which considers the contract term between the movie exhibitor and the 
distributor the revenue corresponding to movie allocation to each week-screen slot 
can be calculated. We propose a Genetic Algorithm based heuristic solution to this 
problem. We compare the improvement in the objective function value (here, total 
revenue) shown by the proposed heuristic over some other heuristics under different 
settings of input factors such as the capacity of screens, the contract term between 
exhibitor and distributor and the number of rapidly decaying movies in a season. We 
find that the proposed GA based heuristic shows an improvement ranging from 3% to 
60% under certain settings of contract terms, type of movies released in the season 
and the capacities of screens. This analysis suggests some actions for the movie 
retailers and distributors that we discuss in managerial implications of the proposed 
research. 
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CHAPTER 2 


LITERATURE REVIEW 

Previous researchers have recognized the inlerl'ace between marketing and operations 
as an important research domain in the quantitative modeling research stream. In the 
operations management literature, Acquilano and (’base (1991, p. 17) mention that 
marketing specialists need an understanding of what the factory can do relative to 
meeting customer due dates, product customization, and new product innovation. In 
service industries, marketing and production often take place simultaneously, so a 
natural mutuality of interest should arise between marketing and OM [operations 
management]. Karmarkar (1996) stresses the need to do more integrative research 
between marketing and operations management. This thesis addresses such needs and 
uses both marketing and operations management tools in solving management 
problems. 

2.1 Shelf Space Management Problems 

Shelf (display) space management is a critical issue to the retailers. First, consumers 
choose from the products that are displayed on the shelf. In this sense, shelf space 
management determines the ultimate profitability of the retailer. Second, most 
retailers have limited shelf space to display their products. Therefore, the choice of 
which products to display becomes an important retailing decision. The decision is 
particularly complex in those industries in which many new products are continually 
introduced in market. Finally, anticipating and adapting to dynamic changes in 
consumer's tastes and demand are key concerns to most retailers. Moreover, such 
decisions have to be made for several product categories. Though previous 
researchers (Corstjens and Doyle 1983; Bultez and Naert 1988; Borin, Farris, and 
Freeland 1994) have addressed some of these issues, they have generally been in the 
context of packaged goods in a supermarket chain setting. However, other industry 
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settings, such as the motion picture (or movie) industry considered in this thesis, pose 
different challenges and intriguing problems. 


2.2 Perishability 

One of the interesting and important attributes of the motion pictures is their 
perishability. In recent years considerable work has been done on the treatment of 
perishability in inventory eontrol (c.g., Abad 1996, .lain and Silver 1994). In inventory 
control, perishability rei'ers to the physical deterioration of units of a product. For the 
movies, we analogously define perishability as the decrease in the value/appeal of a 
movie with the passage of time. As mentioned earlier, previous research in shelf space 
management and dynamic decision making has mainly focused on non-perishable 
goods. In the context of movies, the complexity introduced by perishability in 
dynamic shelf space management problem can be summarized as: which motion 
pictures to choose to show each week and for how long to play them. This decision- 
making problem is further complicated by some other factors specific to the movie 
context. We discuss such factors in detail in the later chapters. 

2.3 Dynamic Decision Making Problems and Markov Decision 
Process 

The dynamic decision making problems fall under a class of problems known as 
sequential decision model (I’ulerman 1994). In these problems, at a specified point in 
time, a decision-maker, or controller, observes the state of a system. Based on this 
state, the decision-maker chooses an action. The action choice produces two results: 
the decision-maker receives an immediate reward, and the system evolves to a new 
state at a subsequent point in time according to an equation of motion determined by 
the action choice. At this subsequent point in time, the decision-maker faces a similar 
problem, but now the system may be in a different state and there may be a different 
set of actions to choose from. The objective of the decision-maker is to choose a 
sequence of actions, which causes the system to perform optimally with respect to 
some predetermined performance criterion. 
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When the decision is assumed to be made at discrete points in time and the 
equation of motion is assumed to be deterministic, the problem is usually modeled as 
a (deterministic) dynamic programming problem (Bellman 1 957). The idea behind the 
dynamic programming approach is to breakdown a large problem into several smaller 
stages -- called separability -- to devise a sequence of simpler optimal solutions that 
lead to the overall optimal solution. It is also possible to write an equivalent linear 
program for a finite state dctenninistic dynamic programming problem (Bertsekas 
1987; Ahuja, Magnanti, and Orlin 1993). This equivalence between dynamic and 
linear programming allows for the problems involving binary decision variables to be 
formulated as integer programs. 

The Markov decision process is one of the sequential decision models 
(Putennan 1994). The MDP considers stochasticity in the transition of system from 
one state to another. Unlike deterministic dynamic programming, the final state of 
system is not detemiined exactly by the action taken by the decision maker. With an 
initial state and action by decision maker, the system has a probability distribution for 
its evolution to the final states. 

2.4 Integer Programming Approach 

A linear programming problem in which some or all the variables must take non- 
negative integer values is referred to as integer linear programming problem. When 
all the variables are constrained to be integer, it is called a pure integer-programming 
problem, and in case only some of the variables are restricted to have integer values, 
the problem is said to be a mixed integer-programming problem. In some situations 
each variable can take on the values of either zero or one; such problems are referred 
to as zero-one programming problems. 

Integer programming is a valuable tool in operations research having a good 
potential for applications. Such problems occur quite frequently in business and 
industry. All the assignment and transportation problems are integer-programming 
problems. In these types of problems the decision variables are either zero or one. 

In all such situations, the decision variable Xj, 

Xj = 1 if j*'’ activity is performed. 
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0 if activity is not performed 

In addition, all allocation problems involving the allocation of men and 
machine give rise to integer programming problems, since such commodities can be 
assigned in integers and not in fractions. 

I'he integer programming problems are solved by either the Cutting Plane or 
Branch and Bound method. In Cutting Plane method by Gomory (1958), we first 
solve the integer programming problem as ordinary L.P. problem and then introduce 
additional constraints one after the other to eliminate certain parts of the solution 
space until an integral solution is obtained. 

In Branch and Bound method, the problem is first solved as a continuous L.P. 
problem ignoring the integrality condition. If in the optimal solution some variable, 
say Xj is not an integer, then 
Xj * < Xj < Xj * +1 

where Xj and Xj * + 1 are consecutive non-negative integers. It follows that any 
feasible integer value of Xj must satisfy one of the two conditions, namely 


Xj<Xj‘ orXj>Xj*+l. 

These two conditions are mutually exclusive and when applied separately to the 
contiiuious L.P. problems, form two dilTcrciil sub-problems. Thus the original 
problem is "Branched" into two sub-problems. Each of these sub-problems is then 
solved separately as a linear program, using the same objective function of the 
original problem. If any sub-problem yields an optimal integer solution, it is not 
further branched. However if it yields a non-integer solution, it is further branched 
into two sub-problems. This branching process is continued until each problem 
terminates with either integer valued optimal solution or there is evidence that it can 
not yield a better one. 

2.5 Parallel Machine Scheduling 

A number of machines in parallel are a setting that is important from both the 
theoretical and practical points of view. From the theoretical viewpoint, it is a 
generalization of the single machine. From the practical point of view, it is important 
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because the occurrence of the resources in parallel is common in the real world. One 
may consider scheduling parallel machines as a two step process. First, one has to 
determine which jobs are to be allocated to which machines; second one has to 
determine the sequence of the jobs allocated to each machine (Baker 1993, Pinedo 
1995). 

In the parallel machine scheduling, the objectives are specified in terms of two 
perfonnance measures, which are as follows. 

2.5.1 Minimizing Makespan 

The makespan Cmax, defined as max (Ci, €„) where Ci,... Cn are the completion 
time of the jobs, is equivalent to the completion time of the last job to leave the 
system. A minimum makespan usually implies a high utilization of the machine(s). 

McNaughton (1959) presented an elementary result for the makespan problem 
when the jobs are independent and preemption is permitted. With preemption 
allowed, the processing of a job may be interrupted and the remaining processing can 
be completed subsequently, perhaps on a different machine. 

The minimum makespan M* is given by 

M* = Max 

m ' 

Where m = number of machines 
n = number of jobs 
tj = processing time of job J 

2.5.2 Minimizing Mean Flow Time 

The other performance measure is minimizing the mean flow time. The flow time is 
the amount of time job / spends in the system. Flow time is given by, 

Fj = Cj - t'l and mean flow time is given by 



where Q is the completion time and rj is the release time. 
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Flow time measures the response of the system to individual demands for 
service and represents the interv-al a job waits between its arrival and its departure. In 
the current study, our objective is neither of the two objectives discussed above, but to 
maximize the total cumulative revenues generated over the time horizon. 

2.6 Decision Support Models in Marketing 

Ovci the last three decades, decision support modeling has flourished in marketing. 
Today one can name many successfully implemented decision support models in 
marketing, such as PROMO'flONSCAN (Abraham and Lodish 1993), Rangaswamy, 
Sinha and Zoltners's (1990) model on sales force restmeturing, SHARP (Bultez and 
Naert 1988), ARTS PLAN (Weinberg and Shachmut 1972), BRANDAID (Little 
1972), CALLPLAN (Lodish 1971), and so on*. Ihese models have addressed various 
aspects in marketing at different levels of a supply chain. For example, CALLPLAN 
efficiently allocates salesforce to customers and products. ARTS PLAN helps 
manager plan a series of performing arts presentations. PROMOTIONSCAN helps 
manager in developing and evaluating short-temi retail promotions. SHARP model 
helps retailers decide on shelf-space allocations, 'fhc SilverScreener model (Swami, 
Eliashberg, and Weinberg 1999), presented in Chapter 4, is similar to 
PROMOTIONSCAN and SHARP models, its aim is to help retailers (of the motion 
picture industry) improve their decision-making. 

Marketing decision support systems (Ml)SS) (Little 1979, p. 9) are intended to 
assist decision-makers in taking advantage of available information. Decision-makers 
should benefit from the availability of more or better data by incorporating the 
infomiation derived from these data into their decision processes (Blattberg and Hoch, 
1990). For this purpose MDSS contain marketing models that make it possible to 
perform so-called what-if analysis (Wierenga et al. 1994). Blattberg and Hoch (1990) 
report that a combination of model and manager often outperforms cither of these two 
alone. Hoch (1994) attributes this to the relative strengths of models, which 
compensates for the relative weaknesses of managers. In general, some of the ways in 
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Many commercial models based on the similar concepts are now also available in the market. 
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which the use of decision support models can aid a marketing executive are the 
following (Montgomery- and Weinberg 1973): 


• Helps to better utili-/e a manager’s judgment, 

• Requires an explicit listing of input assumptions w’hich leads to more informed 
discussion, 

• Provides a method for quick and convenient evaluation of the consequences of 
alternative plans, 

• Allows the emergence of unexpected solutions which open up new areas of 
problem-solving, 

• Expands the range of questions which can be answered by use of the notion of 
derived j udgment, 

• Distills from available data relevant information as in new product forecasting, 

• Provides a basis for relating marketing inputs to market results and, hence, serves 
as basis for marketing planning, and 

• Diagnoses, based on early data, the adequacy of a market plan and locates areas 
needing improvement. 

Inspite of these benefits, one must recognize (he following important points about 
decision support models. First, models are an aid to the decision-maker, not a 
replacement. Marketing models often can help the manager make a better decision, 
but models do not make executive decisions by them.sclves. Second, models should be 
proposed as useful tools to the end-user (i.e., manager), not as academic curiosities. 
Finally, models can be useful to managers in many different ways, and may offer 
opportunities for efficiency in a broad range of managerial activities. 

2.7 Genetic Algorithms 

The Genetic Algorithm approach was first proposed by Holland (1975). The Genetic 
Algorithm seeks to mimic the approach of nature in the evolution of species, that is, it 
is based on the principles of population genetics to guide the search which results in 
“survival of the fittest’. This approach requires the specification of the candidate 


14 



solutions in a binary or a nonbinary string format. These strings are analogous to 
chromosomes in nature. A chromosome in turn is composed of genes, each of which 
can lake a minibcr of values called alleles, rims, for instance, wc have the hair color 
gene at a specific position or locus in the string, which can take a specific allele value 
black, brown or red. In each generation (iteration) the selection of candidate solutions 
to participate in the creation of offspring (new candidates) is based on their ability to 
sur\i\e in the competiti\e environment (that is, its fitness). 'I’he GA approach has 
defined operators to perfomi activities corresponding to the creation of new 
chromosomes as a result of mating and chromosome modification as a result of 
mutation. 

GAs have been applied in variety of fields. Goldberg (1989) has used GA for 
optimization of pipeline systems; Cleveland and Smith (1989) for scheduling flow 
shop releases; Liepins and Potter (1991) use GA for multiple fault diagnosis; 
Balakrishnan and Jacob (1996) use GA for product design; Deb, Chakraborty and De 
(1999) use GA for model based object recognition; Joshi (2000) uses GA in search 
space of hyporcubes. 
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CHAPTER 3 

MOVIE REPLACEMENT PROBLEM: 
MDP APPAROACH 


3.1 Introduction 

Movie making is a large, expensive and risky business especially in North America. 
Hollywood’s major studios produce over 250 feature films each year. These films are 
exhibited on more than 37,000 theater screens in the U.S. and Canada. The average 
movie costs $40 million to produce and an additional $20 million to market. While 
some movies do become “blockbusters,” more than a half of the movies produced by 
Hollywood do not recoup their investment even after their release to foreign markets, 
cable, and television (Vogel 1994). U.S. and Canadian box office revenues have 
continuously been rising in the past decade to exceed $7.7 billion in 2000; however, 
costs have been rising faster than revenues. Thus, improvements in operations and 
decision-making can result in major savings for an exhibitor. 

Every week, motion picture exhibitors have to decide whether to replace the 
movies playing at the screens in their theaters. This is particularly difficult in the two 
major seasons. Summer and Christmas, when a number of movies are released every 
week by the studios-'distributors. Movies are “perishable” products; the box office 
appeal of major Hollywood movies lasts only a short time. This perishability, and the 
complicated contracts between exhibitors and distributors add to the complexity of 
the resulting dynamic shelf-space management problem. We use the term decay to 
explain the intrinsic weekly decline in the box office attraction of a movie playing at 
a theater (Krider and Weinberg 1998). We use another term, aging, to represent the 
decline in the box-office attraction of a movie if there is a delay (e.g., by a week) in 
showing that movie. Aging, therefore, gives rise to an opportunity cost for not playing 
a movie that has already been released in the market. 

The weekly availability of movies and a specialized contract between a 
distributor and exhibitor imply that the decisions an exhibitor takes at one point affect 
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immediate revenue, future decision opportunities, and thus long-term profits. In 
signing a contract to play a film in its theaters, the exhibitor commits to playing the 
film for a minimum obligation period, even if consumer demand is weak. Also, the 
contract splits the box-office gross revenues in a way that favors the distributors in the 
first few weeks of exhibition, but shifts to the exhibitor’s favor later on.' 

Movies decay and age in an uncertain fashion. This adds to the uncertainty in 
demand for movies, and makes the task of their replacement highly challenging. For 
example, if an exhibitor replaces an existing movie with a new one, he/she may be 
worse off than before if the new movie decays very rapidly. In addition, the opemg 
strength of the movie may vary. Finally, distributors release their movies at differ^ 
times during a season; all movies are not available for showing at the same time.' 
However, distributors usually indicate release date of their movie well in advance. 
Thus, it might be beneficial for an exhibitor to continue playing a movie that is not 
doing well, to retain space for an anticipated blockbuster movie to be released stortly. 
In summary, the exhibitor’s movie replacement problem we address is: which motion 
picture to show each week, when to replace it, and if replaced, with what film. As 
discussed below, our primary goal in this chapter is to understand the nature of the 
optimal replacement policies for the perishable products, such as movies, in unc^tain 
and dynamic environments. Accordingly, we focus on the case of a single exhS)itw 
managing a single screen theater and address multiple screen effects in next chapter. 

Vo capture the inter-temporal dependencies of the movie replacement 
problem, we model it as a Markov Decision Process (MDP). At each decision pointy 
in an MDP, the decision-maker observes the system’s state and chooses a cour^ of 


’ For example, the modal contract for one exhibitor across more than 100 movies was for the distributor 
to receive 60%,50%, 40%, and 35% respectively of the box office receipts in each of the first weeks of 
a movie's run. How'ever, if the box office surged beyond a certain prc>spccified level, then the 
distributor would receive 90% of all receipts above that level. In any case, the exhibitor retains all 
concession revenues. For a more detailed explanation of the contract, see Swami, Eliashberg, and 
Weinberg (1999). 

* The multiple screen theater case is a simple e.xtcnsion of single-screen theater case as far as movie 
replacement is concerned. This is because as long as equal screen capacity assumption is made, the 
decision in a multiple-screen case is a vector of movies. This only increases the dimensionality of the 
problem without changing the nature of the model proposed. In different-capacity multiple screen case, 
a modification is required in the definition of revenue generated by a movie on a screen (i.c., minimum 
of potential revenue and screen capacity). This complicates the problem. 1'his problem with the added 
constraint is addressed in next chapter. 
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action. This choice affects future system evolution and reward. An attractive feature 
of the MDP model is that it provides the decision-maker with an optimal policy or 
strategy, that is, a prescription for choosing the optimal action in any possible future 
state. 1 lowever, this feature is also somewhat limiting because it imposes a limit on 
the dimension of the problems that MDP can solve. 

The objectives of this chapter are the following. First, we provide an MDP 
model for the exhibitor’s decision problem. 1 he MDP approach helps us define the 
decision maker as a “smart” manager as one who takes into account uncertainty in the 
environment and adopts a dynamic and long-term view in his/her decisions. Second, 
we investigate optimal replacement policies under different simulated scenarios and 
investigate their managerial implications. Next, we characterize the form of the 
optimal policies that emerge. In particular, we explore the possibility that the optimal 
policy possesses a simple structure in the movie replacement problem. The optimal 
policy results are used as a benchmark for comparing the profitability of simple 
control limit policies and heuristic rules. Finally, we present an example of the real- 
world implementation of the MDP approach. To best of our knowledge, this is the 
first attempt to apply dynamic decision-making techniques in movie theatre 
management. 

The chapter is organized as follows. Section 3.2 presents the modeling 
framework. We present the policies that emerge in various scenarios analyses in 
Section 3.3. In Section 3.4, we compare the MDP model with simple decision 
heuristics. Section 3.5 presents an example of the application of the MDP model 
approach in a real-world setting. We conclude in Section 3.6 with a discussion of 
limitations and extensions. 

3.2 The Markov Decision Process Model 

3.2.1 Movie Replacement Problem 

Every week, motion picture exhibitors make decisions regarding tlic replacement of 
the movie playing at their theater. From a modeling perspective, the key features of 
the movie replacement problem are decay and aging of movies, uncertainty of 
demand, the specialized contract and the fact that all movies are not available for 
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showing at the same time. 'I'he movie replacement problem is somewhat analogous to 
the equipment replacement problem (Derman, 1963) with the current exhibited movie 
occupying the role of the deteriorating machine. 

To model the deterioration of a movie at a screen, we introduce a parameter, 
the rank, of a movie. We define the rank as its relative position with respect to other 
currently available movies in terms of box-office gross revenue. We assume that in 

any week, a movie can be in one of the ranks 1 c, where 1 denotes the highest rank 

and r, the lowest’’. We further assume that knowledge of the rank of a movie provides 
complete information about its expected gross revenue. In other words, different 
movies with the same rank have the same expected gross revenue. Previous 
researchers have also used the notion of ranks to characterize performance of a movie 
(De Vany and Walls 1997). 

We now describe the timing of events in the exhibitor’s decision problem (see 
Figure 3.1). Each Monday morning, after observing the weekend demand, the 
exhibitor must decide whether to replace the current movie (say, a) and, if so, with 
which available film (say, h). In making this decision, the exhibitor must take into 
account the potential revenues of the cuiTcnt and replacement movies, the contract, 
and the release dates of replacement and future movies. At the beginning of the 
season, the release dates of all movies for that season are usually known. Randomness 
enters the decision problem through the reward streams (attendance) of the movies, 
which follow a probabilistic pattern. We emphasize that if a movie has been released 
and not shown by the exhibitor, its appeal may decline. 
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Figure 3.1: Timing of Events in the Movie Replacement Problem 


Start of Decision Start of 

Week i Epoch/ Wcckf/+/j 

Weekend demand 


FRIDAY- MONDAY- FRIDAY - 

Movie ii Movie a playing - Movie h playing 

playing Choose movie h to play 

on the coming Friday 

We assume that the exhibitor has sufficient information to determine a 
movie’s rank (an indicator of its expected attendance) for the upcoming decision 
cpooli. This appears to be a reasonable assumption. Swami, liliashberg and Weinberg 
(1999) report that week ahead forecasts were correlated at 0.96 with actual box office 
revenues for the movie theater analyzed in their paper. 

In the case of an already released movie, the exhibitor might know its rank 
based on its attendance pattern at this theater in previous weeks, because an additional 
print might be playing at another theater or from information about its nationwide 
performance. In case of a new movie, sources of information for the exhibitor to be 
confident of a movie’s opening include advertising levels, word-of-mouth, private 
market research, and so Ibrth. 

The transition between ranks over time is represented by a Markov process, 
which is designed to replicate the empirical patterns found in Jedidi, Krider, and 
Weinberg (1998). Also, based on Jedidi, Krider, and Weinberg’s (1998) results, we 
assume that after the release of the movie, its rank is non-improving, that is, it does 
not make transition to a higher rank 

We now formally state our modeling assumptions: 

1) The theater has one screen. 

2) The opening rank of a movie is known before its release. 

3) Replacement decisions are made on a weekly basis. Further, the replacement 
decision for the coming weekend is made on the previous Monday, and delivery 
of the replacement movie occurs instantaneously. 

Note that we augment the rank designation when wc formulate our model below. 


20 



4) The release dates of the movies considered during the planning horizon are 
deterministic and known in advance. 

5) Movies once replaced at a theater arc not available for subsequent .screening at the 
same theater. 

6) The probability of a transition of a movie between ranks is stationary over time 
and independent of the other movies. 

Assumption 3, 4, and 5 are consistent with industry practice. Assumptions 1 and 6 
have been made for the tractability of the analysis. Assumption 6 might be limiting 
because of the following reasons. The independence between the transition 
probabilities of movies may be questionable in some cases, because decay of a movie 
may affect that of another. Second, the transition of the movie between the ranks may 
vary depending on whether the movie is in the early or later part of its run. We make 
this a.ssumption to keep the analysis simple at this level of generality. The relaxation 
of the assumptions will provide attractive future research opportunities. 

3.2.2 Model Formulation 

A Markov decision process (Puterman, 1994) often referred to as an MDP, is a model 
for a stochastic sequential decision process. To formulate an MDP we specify 
decision epochs, states, actions, rewards and transition probabilities. In our 
application, the exhibitor notes the movie currently playing, its length of play and the 
ranks of all movies that could be chosen to replace it. The exhibitor then decides 
whether or not to replace the movie and if so, with which available film. As a 
consequence of this decision the exhibitor receives a random stream of daily revenues 
throughout the current week and all movies, which have already opened age and 
decay in a probabilistic fashion. By solving the MDP, the exhibitor obtains a policy 
or contingency plan which specifies which action to choose each week given any 
possible configuration of movie being played, its length of play and ranks of available 
replacements, so as to maximize expected total revenue over the planning horizon. 

Before fonnally defining the model, we introduce the following notation. 
Recognizing the complexity of the notation, we illustrate and describe it following the 
definitions. 



W — length of planning horizon, in weeks, say summer season, 

M — set of movies available during the planning horizon, M={mi,m2„in3 hin}, 

N — total number of movies considered during the planning horizon, A^=|M|, 

R — set of ranks of a movie, ^={ 0. 1 2. r +1 

— set of all possible combinations of ranks of the N movies in the set M, 
r - a typical vector of ranks of all movies 
OPD — the length of the obligation period,'^ 

?}i — index number of the movie playing at a decision epoch, in e M, 
u — total number of weeks the movie in has played before a decision epoch, 

Z - vector of ranks, at a decision epoch, of the movies considered during the planning 
horizon, 

£'[A^i?( 5 .a)] -- expected net revenue to the exhibitor at a decision epoch given state s 
and action a, 

Pj(h\h) — probability of transition of movie j, when it is available, from the rank ij to 
ij), where ii.i2 e {l,...,z}, 

P'j (0 -- the probability that movie j opens in rank i ,i e{ 1 . .,z}, 

tm, h’ " release dates of the movies in A/ during the planning horizon. 

We now explain the notation in detail. The length of the planning horizon, W, 
is finite and measured in weeks.'*’ An example of a planning horizon is the summer 
season, which extends in North America from the last week of May (Memorial Day in 
the U.S.) to the first week of September (Labor Day in the U.S.). The set of movies, 
M, consists of a list of the movies available at any time throughout the planning 
horizon. For example, if three movies, a, b, and c, are available, then the first movie, 
nil = a. Similarly, 1112 - b and 1113 = c. Consequently, M = {a,b,c) and the total 

■* The length of the obligation period is proposed as movie-independent for the case of simulation 
analyses to be introduced later. However, the model can be readily extended to make this parameter 
movic-dcpcndcnt. 

^ We assume that the last decision is made at decision epoch W. Since we do not know future film 
availability beyond the planning horizon, we choose to fomiulate a finite horizon problem. The 
consequence of doing this is that the terminal condition might affect optimal policies in some cases. 



number of movies, iV, equals three, l^he set of movie ranks, R, lists the possible ranks 
an available movie can assume from 1, the best rank to r, the worst rank. 1^'or example, 
in a ranking system {High, Medium, Low}, High ~ 1, Medium ~ 2, Low = 3, and 
consequently, z = 3. We augment the ranks, 1 to z, by two levels, 0 and z+1, to denote 
the movies that are not available for play. If a movie has not been released, its rank is 
0, and if it has been replaced after having been shown by the exhibitor, then its rank is 
set to z+1. The set, , gives all possible combinations of ranks of N movies. For 
example, if iV = 3 and i? = {0,1, 2, 3}, then = {(0,0,0); (0,0,1); (0,0,2); (0,0,3); 
(0,1,0); (0,1,1);...; (1,0,0); (1,0,1);...; (2,0,0); (2,0,1);...; (3,0,1 );...;(3,3,3)}. The 
vector of ranks, 2, of various movies has dimension N and denotes the ranks of 
various movies at a decision epoch.^ It is clear that Z can assume any value in at 
any decision epoch. In the above example, if Movies a and b have not been released, 
and Movie c is available in Rank 3, then 2 = (0,0,3) e R^. 

We now provide our MDP formulation of the movie replacement problem. We 
explain each component of this apparently complex model by an illustrative example 
in Section 3.2.3. The approach to evaluating a policy for implementation is discussed 
in Section 3.3. 

Decision Epochs (beginning of every week): 

r = {i,2 «' jk},w<co. 

States (the movie playing, its play length after obligation period, and the ranks of all 
the movies at a decision epoch): 

5= {/n,«,2}, m e M, n s {-OPD, -OPD+1,..., 0,1,2,..., W-OPD},ZeR^. 
Actions (continue or replace a movie when possible): 

f 111 if n < 0, or 2^ e {0,z + 1} ybr a//y e M - {/«} 

. 2 ^ g ^Q - + ij' „ > 0 aiul Z ^ € {0,z + 1} for all j & M - {m} 


However, we did not observe any systematic effects by varying the planning horizon length in this 
study. 



Rewards (expected net revenue to exhibitor in week vv): 


rH(s,a) = £lNR(s.a)], w = s e S, a € As . 

/•(('+/ (s) = 0, s e S. 


Transition Probabilities: 

p„.(s: 1 s/. a) = 

'n (Z IZ t if {a = m,=m^,n, = n, + \,Zj.>Z-^>0,foraUjeM) or 

,.2l y.i)> ^{a = nu^m^,n,=-OPD,Z„^^, = z + lZ..>Z.,>0,foraUjsM) 
I 0, otherwise. 


where 


pAZj, 


= 0|Z 


J.l 



if M' < r - 1 
otherwise, and 


(movie not released) (3.1) 


=A-|Z,,, =0) = /;;"(A'), u- = /,-l, A-e|l,...,;j, Z, ^0 V u->/, -1, 


(movie to he released next week) 
Pj{Z ^ 2 ~ M 2^ ^ + 1 ) = 1 (movie already played and replaced) 


lor 


(3.2) 


si, s: 0 S,si - (mo ///. //), - (m:, Zj); a e /l.v; u' - U..., IV. 


3.2.3 An Illustrative Example 

We begin with the hypothetical movie release scenario shown in Figure 3.2. We 
assume three movies with two possible ranks for each. 


* Zjj denotes the rank of movie J at decision epoch t. 



Figure 3.2: Movie Release Pattern for Illustrative Problem 


Movie 

Released a 


Week 1 2 3 4 5 6 7 8 



Figure 3.2 shows that Movie a is released and available in the first week. Movies c 
and b become available in Weeks 4 and 5 respectively, that is, t„=l, t/,=5, and tc=4. 
We choose the following values for the other parameters; W = M = {a,b,c}, N =2>, 
/?= {0,1, 2,3}, OPD = 2, 


a).3 0.7'' ^ _ ^0.5 


0.5' fO.l 
P = 

1 / ^ V 0 


o.W'- 

, and 

1 > 


p': = (0.2 0.8), pt = (0.7 0.3), p:: = (0.8 0.2) . 


In this example, the length of planning horizon is eight weeks. The obligation 
period is 2 weeks for all movies. The movies can be in the ranks 1 (high) or 2 (low) 
when they are available. Rank 0 implies that a movie has not been released and Rank 
3 implies that it has been replaced after having played. The transition and the initial 
probabilities are specified for ranks 1 or 2. We express them in matrix format for 
simplicity. For example, the matrix P„ has components PnOV) which denote the 
probability that movie a occupies rank / at the next decision epoch given that it 
occupies rank j at the current decision epoch. For example, Pa(2\l) = 0.7 means that 
there is a probability of 0.7 that movie a declines from rank 1 to rank 2 in any week. 
These probabilities imply that Movie a is the weakest of the three movies. This is 
because its initial probability of opening in Rank 1 , the better rank, is low (0.2), and 
even if it does, there is a high probability (0.9) that it decays to Rank 2. The expected 
net revenue generated by a movie at a decision epoch depends on its rank and the 
number of weeks it has played before that epoch. The expected net revenue to the 
exhibitor, presented below, has been generated by assuming common contract terms, 
which are representative of the industry practice. 



^ Rank’ ' 1 ' ' 2 3 4 5 ' 6’ 7 8 ’ 

_\\Vcck , ^ _ 

1 ’ 290 ' 290 ”290’ ’ 290 ' 290 ' 290 290 290 : 

2_ 60 ; 70 80 90 90 : 90 ’ : 90” ’ 90 

An explanation of the MDP model Ibllows. Suppose that we observe the 
system in Week 6 and find that Movie c has been playing for two weeks after 
replacing Movie a in Week 4. Therefore, m6 = {c}, and «i 5 = -2 + 2 = 0, because the 
obligation period is two weeks. Further, suppose that the ranks of Movies c and b are 
1 and 2 respectively in Week 6, which implies that r* = (3,2,1 ). The rank of Movie a 
is 3, because it has been replaced. Thus, the state variable, S6, = (c,0,(3,2,l)). 

The action specifies whether to play the current movie for an additional week 
or to replace it and if so, with which other available movie. This decision option 
depends on which movies are available for replacement, and whether the movie 
currently playing remains in obligation period. The only action is to continue playing 
the current movie if either there is no movie available for replacement, or the current 
movie is in its obligation period. If both of these conditions are false, then the action 
set consists of both continue and replace decisions. In the example problem, the 
feasible action in the first three weeks is a, because no other movie is available. In 
Weeks 4 and 5, since Movie a is replaced (its rank is 3), and Movie c is still in its 
obligation period {tu and ns < 0), the set of feasible actions is {c}. From Week 6 
onwards, because Movie c has played longer than its obligation period, the set of 
feasible actions is {Ac}. 

The reward earned by the exhibitor at decision epoch vv is the one-period 
expected net revenue received by the exhibitor by taking action a„. in state We 
assume that the rank of a movie determines its expected gross revenue. We set the 
reward at period W+1 (i.e., salvage value) to zero*^. Suppose Movie a is playing in 

’’ For the sake of simplicity, we do not provide the reward function in this paper, but simply note that it 
depends on the contract terms, c.xpcctcd gro.s.s revenue, and c-speeted eoncessioti profits, fhe c.xact 
expression of the reward function is similar to the one used by Swami, Eliashberg and Weinberg 
(1999). When the box office returns are sufficiently high, the exhibitor retains a constant and not an 
increasing share of these receipts. 

* Choice of terminal reward can have a significant effect on the optimal policy. Alternatively we could 
set the terminal reward equal to the revenue from playing the movie at the last decision epoch 
indefinitely. 



Week 5 (i.e., it has played for a total of four weeks) and the ranks of Movies a, b and 
care 1, l,and 2, respectively. Thus, 55 = {a, 2, (1,1,2)}. Then, if either movie a or Ms 
chosen, tlic one-period immediate e.\pected reward to llic c.xhibiU)r is $290; if Movie c 
is chosen, the reward is $60. 

I'he probability transition function, /j„(.v: |.s/, a), specilies the probability that 
the system will be in the state S 2 ={m:. 112 , 22 ) at the next decision epoch given the 
current state of the system is ,v/=(m/.»/.Z/) and action a is taken. The last component 
of the state variable, Z,- contains the extended ranks of all movies; as noted before, 
those that have yet to open have rank zero and those that have already been shown 
have rank 3. Since the transitions of different movies between various ranks are 
independent of each other, the transition probability between rank vectors is simply 
the product of the individual movie transition probabilities. If the action is to continue 
with the current movie, then the movie index (i.e., the first element of the state 
variable) remains the same in the next period. The play length of the movie increases 
by one period. All the movies that have been released either stay in the same rank or 
make transitions to a low'er rank according to their transition probabilities. If instead 
we choose the action to replace the current movie by another available movie, the first 
element of the state variable in the next period becomes the movie index of the new 
movie; the second component, the play length, is set to OPD denoting that the movie 
has entered its obligation period. The rank of the current movie m\ is set to z+1 in the 
next period. The rest of the available movies make transitions according to their 
respective transition probabilities. Some movies may not be available at a particular 
decision epoch because either they have not been released by that week or they have 
been replaced after having played. The probability expressions for these cases are 
given by Equations 1 and 2 respectively. Equation 1 implies that if a movie has not 
been released (i.e., its rank is equal to zero), then its rank remains zero with certainty 
in the next period until a period before its release. At a period before its release, the 
movie has an initial probability of opening in the next week in ranks 1 to z. Equation 2 
implies that once a movie is replaced (i.e., its rank is equal to z+ 1 ), it is never 
available for screening again. Suppose Movie a is playing in Week 6 and all the 
movies are in Rank 1. Therefore, si = {a,3,(l,l,l)}. If we choose action {£>}, then the 


probability of system making transition to .9: = \h, -1, (3,1,2)} is 1 x 0.5 x 0.9 = 0.45. 
Alternatively, suppose Movie a is playing in Week 3 and its rank is 1, then the state is 
{a,0,( 1 ,0,0)J . Movie c becomes available in liie next week. The probability of 
transition to {a,l,(l,0,l )}, therefore, involves the initial probability of Movie c 
opening in Rank 1 and movie a stating in Rank 1, and equals 0.3 x 1 x 0.8 = 0.24. 

3.2.4 Determining Optimal Policies 

By a decision nile, we mean a function that specifics the action to use in each slate at 
a specified decision epoch. We use the notation to represent a decision rule. It 
gives the action chosen in state s at decision epoch vv. A policy is a sequence of 

decision rules, one for each decision epoch. We represent policies by n = (di, d: 

da ) and compare them on the basis of their expected total reward over the planning 
horizon. Let v’‘(s) represent the expected total reward if policy n is used and the 
system is in state s at the first decision epoch. Then, 

v'' is) = E: (£ T„, ) + r,,,, ( ) } (3.3) 

VV=rl 

where Es’' represents expected value under the policy n when the system starts at 
decision epoch 1 in state s and Xw and T„. are random variables, which denote the state 
and action at decision epoch vv, respectively. An optimal policy, n, with the largest 
expected total reward, defines the value, \ '(s), of an MDP model through 

v(i-) = v {s)= {v (i-), .vei’} (3.4) 

TT e n 

where El denotes the set of all policies. In the next section, we compare and evaluate 
various alternative scenarios on the basis of the value of their MDP model. We define 
a smart manager as one who uses the optimal policy recommended by the MDP 
model and achieves the maximal total expect reward. We use the backward induction 
algorithm (Puterman 1994, pp. 80-82) to determine an optimal policy and value 
function. 

To illustrate these concepts and the benefits of this approach, we interpret the 
optimal policy for the example problem shown in Table 3.1. 


Table 3.1: Optimal Policy for the Illustrative Problem in Section 3.2.3 


Decision epoch 

T 

State variable {.v„.) 


Optimal 

(w) 




decision 





rule 


Movie 

playing 

Number of weeks 
played 

Ranks of 
movies 

\ 

1 to 3 

a 

’ Corresponding to ' 

0 or 2,0,0) 

a 

_ 


decision epoch 

— , - 


4 

a 

2 

(2,0,1) 

c 




Any other 

a 




combination 


5 to 8 

a 

Corresponding to 

(2,2,1) 

c 



decision epoch 






(2,T,i or'2y 





Any other 
combination 

a 


b 

Corresponding to 

(3,2,1) 




decision epoch 

Any other 
combination 

b 


c 

Corresponding to 

(3,1,2) 

b 



decision epoch 

Any other 

c 


combination 


The table shows the state variable, 5 ,,., and optimal decision rule, be 

followed in Week w of the planning horizon. 

Weeks 1 to 3 - (Movie a is playing at the beginning of the horizon and no other 
movies are released until Week 4.) The optimal policy chooses the only feasible 
action a. 

Week 4 - (Movie c becomes available.) Replace Movie a by c only if c opens in Rank 
1 and a has decayed to Rank 2 , otherwise play Movie a again. 

Week 5 to 8 - 

a) If Movie a is playing in Rank 2, then replace it by Movie c if both movies a 
and h are in Rank 2, otherwise replace it by Movie h if Movie a has decayed to 
Rank 2 and Movie h is available in Rank 1. In the latter case, the rank of 
Movie c does not affect the optimal policy. This is because even if Movie c is 


on 



available in Rank 1 , it is highly likely (with probability 0.9) to decay to Rank 
2 by the next week. 

b) If Movie h (or c) is playing and beytnid its obligation period, then replace it 
by Movie c (or h) only if the replacement movie is available in a better rank. 
This implies that if the ranks of the two movies are equal, then it is optimal to 
play Movie h (or c) again. 

3.3 Simulation Analysis 

The optimal policy for an exhibitor to follow depends on the quality and quantity of 
movies available and their release timings. We conduct an experimental study to 
investigate the nature and magnitude of these effects on optimal policy. 

3.3.1 Model Parameters 

The input parameters required to generate problems for various simulation analyses 
are as follows: 

(i) length of planning horizon 

(ii) number of ranks 

(iii) quantity of movies (i.e., total number of movies considered) 

(iv) quality of movies 

(v) release dates of movies 

(vi) expected revenue a movie generates in a given rank 

(vii) probabilities of transition of a movie between ranks (including initial 
probability) 

To generate different simulation problems, we assign fixed values to the first two 
parameters, which are common to all problems. We set the length of the planning 
horizon to be 8 weeks (e.g., representing the two peak months of Summer season). 
We consider a three-rank system (e.g.. High, Medium, and Low) for the box-oflicc 
performance of movies. The parameters quantity (iii), quality (iv), and release dates 
(v) are specific to each simulated problem scenario, and are discussed in a later 
section. We now discuss a movie industry based study by Jedidi, Krider, and 


in 


Weinberg [JKW] (1998), which helps in the operationalization of these and the 
remaining parameters (vi. vii) of the problem. 

JKW modeled the revenue pattern of 102 major motion pictures released in 
North America during the period December 1990 to April 1992. They used an 
exponentially decaying function to model the box-office revenue in week t as given 

<X +£ 

by " , where a denotes a movie’s opening strength and /3(<0) its revenue 

decay rate. JKW used cluster analysis to segment the movies on the basis of the 
estimated a and /3 values, and identified four clusters of films which they referred to 
as Type 1: Hollywood Heroes (movies that open strongly but decay quickly, e.g.. 
Awakening, Sleeping with the Enemy), Type 2: Mega Movies (movies that open 
reasonably strongly and decay slowly, e.g.. Home Alone, Pretty Woman. Silence of the 
Lambs), Typ® ^'^st Fades (movies that open weakly and decay very sharply, e.g.. 
Bad Injluence, Lord of the Flies, White Palace), and Type 4; Fair Flicks (movies that 
open weakly and decay sharply, e.g., Joe Versus the Volcano. Internal AJfairs, 
Postcards from the Edge). The parameters a and p for the four types of movies are 
shown below. 

Type 

1 2 _ _3_ , ' 4" 

a ' -1.484 -1.6l'4 ' -1520 ^ '-2.255 

j3 -0.224 -0.110 -0.439 -0.258 

Based on the above study, we now discuss the generation of the other parameters of 
the simulated problems. 

3.3. 1.1 Exhibitor’s Revenue 

In the proposed model, the rank of a movie detemiines expected box-office gross 
revenue. To operationalize the three ranks, 1) High, 2) Medium and 3) Low, we first 
assumed a range of weekly box-office revenue from 0 to in excess of 1000 units. To 
generate the three ranks, we assumed threshold values at 1000 and 500. Thus, a movie 
which generates more than 1000 in revenue is classified in the high rank, a movie 
which generates less than 500 units is classified in the low rank, and the remaining 



cases are classified as medium. JWK’s model is normalized so that the weekly 
revenues in Ranks 1, 2, and 3 are 1250, 750, and 400 units, respectively. 

However, even if two movies are in the same rank, the revenue the exhibitor 
earns will differ according to the contractual sharing terms, which vary by the number 
of weeks since a movie's release. To calculate the exhibitor’s share of the box-office 
revenue generated by any movie in a given rank and week (refer to Table 3.2), we 
chose the sample contract terms used in Swami, Eliashberg, and Weinberg (1999). 

Table 3.2: Expected Weekly Net Revenues to the Exhibitor 


Rank Run-length 



1 

" 2 

' T- 

3 

4’ 

5 

6 ’ 

7 

8 

i 

■ ‘440 

"440" 

'""440 " 

440 

" ' 4^ 

440 

440 ' 

, 440 

2 

300 

360 

390 

390 

390 

390 

390 

, 390 

3 

160 

180 

220 

240 

240 

240 

240" 

'"~240 


3.3.1.2 Initial and Transition Probabilities 

Initial and transition probabilities are discrete time representations of the exponential 
declining model. To estimate the values of initial and transition probabilities of each 
type of movie, we use their respective a and /3 values and repeatedly generate weekly 
revenues by varying t from 1 to 50 in the exponential demand model. To induce 
randomness, we use the error estimates given in Sawhney and Eliashberg (1996)’s 
study on forecasting of movie revenues. 

The initial probability of a particular movie type to open in a particular rank is 
operationalized by the proportion of time the movies (for a large generation of 
movies) of this type open in that rank in the first week (refer to Table 3.3 (a)). 
Similarly the probability of transition of a movie, from rank, say r/, to another rank, 
say is operationalized as the proportion of times the movie decays from rank r/ 
from rank rj among all transitions from rank r/ (refer to Table 3.3 (b)). 

For example, JKW’s clustering results show that Type 1 movies open stronger, 
but decay quicker than Type II movies. This pattern is preserved in the probabilities 
generated since Type I movies never opened in Ranks 2 or 3 (corresponding 



probability = 0), while 5% of the Type II movies opened in Rank 2 (probability = 
0.05). On the other hand, the probability of a l ype 11 movie retaining ifsclf in Rank 1 
is higher (0.8) than that of a Type I movie (0.7). 

Table 3.3: Probability Matrices for Scenario Analyses 

(a) Initial Probabilities ' 

Type of movie Rank 

1 2 ■ 3 

I ’ 1 ' 0 ‘ 0 

-,.U .. 0.05 ' ' o ' 

ill "O ' 0.56 0.44 

IV ' 0.07 0.93" 0 


(b) Transition Probabiiijlies_ 

Type I movies 

Rank 1 2 3 

“ 1 0.7 ' 0.3 0 

2 0 ' 0.7 ^ 0.3 

3 _ _ _ ' y 0 „ _ 0 _ J ; _ \ 

Type il moyi^ _ 

RaiTk ' f" “ " 2 ' '3 

1 ' 0.8 ' 0.2 ~ 0 
2 ' 0 0.8 ' 0.2 

__ 3 ' ’() / _ 1 . 1 

Type III movies 

Rank ' 1 2 3 

1 0 o' 1 

2 0 0 ’ 1 

' 3 ' ' 0 ■ ' 0 ^ r 

Type IV movies 

_Rank ^ 1_ " i ""' 3 

T ' “b J ~ Ib., _ o' 

2 0 ' 0.6 ' 0.4 

' "3 ' '■ o' ""' O ■ \ 



3.3.2 Simulation Setup 

We study the ettect of quantity (higli/low), quality (high/low), and release timings 
(high/low early release of quality movies) on the exhibitor’s expected revenue using 
in a complete 3 (treatment) X 2 (level) factorial design. 

3.3.2. 1 Design of experiments 

Wo examine two levels, high and low, ol the quantity elfect. In the low quantity case, 
we assume that a single movie is released in every week of the planning horizon (i.e., 
an 8-movie problem). In the high quantity case, we assume that two movies are 
released in every week of the planning horizon (i.e., a 16-movie problem). 

The quality of the movies is operationalized as the number of high quality 
movies (Type 1 and 2) in a given scenario. To randomize the assignment of a movie 
to a movie type, we use the relative proportions of Type I (19%), II (7%), III (38%), 
and IV (36%) movies in Jodidi, Kridcr, and Weinberg’s (1998) sample. Specillcally, a 
random number is drawn among integers between 1 and 100 for each movie. Then 
the movie is assumed to be Type I if the random number is between 1 and 1 9, Type II 
if it’s between 20 and 26, Type III if it’s between 27 and 64, and Type IV if it is 
between 65 and 100. 

We found in a preliminary analysis that when there are three or more 
high quality movies (Types 1 and 2) in a scenario, their effect on the value is much 
more pronounced than the rest of the movies, 'fherefore, a scenario in which we find 
three or more high quality movies is classified as high quality scenario, otherwise it is 
classified as a low quality scenario. 

To operationalize the release date effects, we define an early release factor 

(e.r.f.), which indicates the proportion of high quality movies that are released in first 

four weeks (i.e., first half of the season). 

„ , „ , Numberof high quality movies released in Weeks 1 to 4 

Itarly Release factor = — — r ^ — : -. — 

I'otal number ol high quality movies in the planning horizon 

Using a cut-off value of 0.5, we classify a scenario with early release factor of greater 
than 0.5 as a high early release effect scenario, otherwise as a low early release effect 



scenario. We use a complete 3-way 2-level factorial design for our experiment. That 
is all 8 combinations of the three factors are evaluated. A summary of the resulting 
designs is provided later in Table 3.4 (a). 

3.3.3 Methodology 

Using the above scheme of simulating problems, a total of 66 problems were 
generated. Based on the input data, the problems were assigned to their respective 
type of experiment. The MDP model of Section 3.2.2 was coded in C language and 
run on an Intel Pentium III class computer. The value evaluations were done 
according to the cumulative expected value criterion presented in Section 3.2.4. 
Although a small-size problem (e.g. 3 movies, 3 ranks, and 8 weeks) would take only 
a few seconds to run for the MDP algorithm, the time taken to solve a problem 
increases drastically with problem size, especially with the number of movies and 
ranks. In the current analysis, time taken to solve the MDP problems ranged from 3-5 
minutes for low quantity cases to 40-45 minutes for the high quantity cases.’ 

Due to the limitation concerning dimensionality in MDP problems, we 
adopted several state-reduction mechanisms, to speed the execution of the algorithm 
without impacting its logic. For example, it was observed that the high quality movies 
(Types I and II) dominate the value of the program. Therefore, we reduced the 
number of movies in the consideration by including low quality movies (Types III and 
IV) only if they have not been older than four weeks (i.e., half the season). Similar 
other mild assumptions were used after checking the robustness of the results. 


’ The time taken to solve a problem depends on the complexity of the problem determined primarily by 
the number of movies and ranks. For example, in a typical low quantity case (8 movies, 3 ranks, 8 
weeks), the number of movie-rank combinations to bo considered for each possibility of the state 
variable are 3®, or 6561. Since these combination are for each of the eight possible movies playing and 
their run-lengths (an average of four weeks), the total number of records examined in a typical stage of 



3.3.4 Simulation Results 


Table 3.4 summarizes the simulation results in several ways. Table 3.4(a) gives the 
mean (and standard deviation) under each factor level combination and Table 3.4(b) 
gives the mean and standard deviation for each main clTect. fhe ANOVA results in 
Table 3.4(c) show that the interactions are not significant and that all three main 
eflects (Quantity, Quality and harly Release Factor) are significant at the A)5 level. 


Table 3.4 (a): Mean Expected Cumulative Revenues for Experimental Setup 


Experiment 

Quality’’ 

Quantity^ 

Early 

release 

factor’’ 

Number of 

problems 

Mean ' 

expected 

cumulative 

revenue* ’ 

1 

1 

0 

0 

6 

1757 (298) 

2 

1 

0 

1 

6 

1939 (419) 

3 

r 

1 , 

0 ’ 

6 

2191 (69) 

4 

1 

1 

1 

6 

2213 (72) 

5 

0 

0 

. 0 

6 

1247(247) 

6 

0 

0 

1 

6 

1901 (170) 

~ 1 

0 

l‘ 

0 

6 

1631 (94) 

8 

”o 

' 1 

1 

6 

1 854 (422) 


Totals 48" V 

a ~ High (1) setting implies greater than 2 high-quality movies in a season, otherwise low (0) 

setting 

b - High (1) setting implies 16-movie scenario, Low (0) setting implies 8-movie scenario 
c - High (1) setting implies an early release factor greater than 0.5, otherwise low (0) setting 

* Standard deviation values are shown in parentheses. 

** The total number of problems is 48 as 1 8 out of the randomly generated 66 problems 
could not be assigned exactly to an experiment type. 


backward induction algorithm arc of the order of 209952, However, a similar number for a typical high 
quantity case is of the order of 2754990144. 



Table 3.4 (b): Mean Expected Cumulative Revenues 
for Different Settings of Parameters 


Factor 

Setting 

Mean Expected 



Cumulative Revenue* 

Quantity 

High 

1972 ( 164 ) 


Low 

171 1 ( 283 ) 

Quality 

“Higii 

2025 ( 214 ) 


Low 

. 1 658 ( 233 ) 

1 

Early release 

High 

1976 ( 270 ) 

factor 

Low 

1706 ( 177 ) 


Standard deviation values are shown in parentheses 



Table 3.4 (c): ANOVA 

Results 


Source of 
Variation 

Sum of 
squares 

Degrees of 
freedom 

Mean square 

Fo 

Quality (A) 

1618030 

1 

1618030 ' 

11.1135 

Quantity (B) 

876680 

1 

876680 

6.02* 

Early release 

820273 

1 

820273 

5.63’ 

factor (C) 





AB 

103253 

1 

103253 

0.709 

AC 

339495 

' 1 

339495 

2.'33 

BC 

' 261528 

" 1 

261528 

1.79 

ABC 

'55284 

“ l' 

55284 

0.379 

Error 

2768239 

40 ' ' 

58898 


Total 

6842782 

47 



• 

Significant at 5 percent 





It can be concluded from the above results that quality of movies released is a 
major factor affecting the value (expected total revenue) the exhibitor achieves in the 
planning horizon, 'fhough tlic on’ccls of quantity and early release lactors are also 
statistically significant, the change in the value is the greatest when the quality 
changes from a low to high setting. The standard deviation is low for the high setting 
of quantity and quality factors. This is because in such cases a few movies, usually of 
high quality, dominate the exhibitor’s policy. However, if such movies are available 
in the early half of the planning horizon (a high early release factor scenario), then the 
deviation in mean response values would be higher depending on when the high 
quality movies are released during Week 1 to 4. Note that even though all interactions 
are insignificant, the Quality-Early Release Factor Interaction has the greatest p- 
value.'” 

To further investigate the effects of these variables on exhibitor’s value, we 
conducted a simple linear regression analysis on the simulation data. In this analysis, 
we decompose the effect of quality into two components: {A) mean opening strength 
of the season (i.e., planning horizon), and (5) mean retention strength of the season. 
The first variable is a measure of how well the movies opened on the box-office in a 
particular season. The second variable is a measure of how well the movies were able 
to retain themselves (i.e., low decay rate) at the box-office. The other two variables 
considered are number of movies {N, quantity effect), and early release factor (/i). 

Mean opening strength of the season is operationalized as the mean of the 

C£ + j3 £ 

opening strengths (a function of a in the exponential model A6 ) of all 

movies that are released in the season. Similarly, mean retention strength of the 
season is operationalized as the mean of the decay rates (a function of p in the 

exponential model € ) of all movies. Specifically, if N movies are released, 

Kridcr and Weinberg (1998) study the effects of opening strength and decay rate on movie release 
strategies; consistent with their findings an empirical regression-based analysis of the data in this 
experiment found that both opening strength and decay rate were significantly related to exhibitor 
value. 



then A and B are quantified as follows. 
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The resulting regression model with exhibitor’s value, V, as the dependent variable, 
and B, N, and E as the independent variables, is given by the following equation. 

V = +b2* A + b^* B + b4* N + b^* E + s {3.5) 

The model results for the simulation data are given below (t-values are in 
parentheses). 

^ = - 3099 + 9376 . 81 *^ + 4257 . 09*5 + 50 . 42 * + 191.77 *5 

(-2.7) (3.38) (2.26) (6.73) (2.14) 

=0.87 

A 

where V is the mean exhibitor’s value for the season. 

The model provides a good fit to the simulation data (R^ = 0.87). The significant 
values of the regression coefficients corroborate the results achieved earlier. In 
absence of a composite measure of quality, the effect of the quantity is the strongest. 
However, the results on the two measures of quality further information regarding 
opening and retention strengths. The results indicate that in the simulated problems, 
the effect of the opening strength is slightly stronger than the retention strength of the 
movies in a season. This could be due to the relatively smaller length of the planning 
horizon, which may also explain the lower (although significant) effect of the early 
release factor. 

3.3.5 Discussion 

Intuition suggests that the exhibitor should be better off in a high quantity than a low 
quantity scenario. The reason is that a high quantity scenario offers the exhibitor 
increased flexibility in scheduling movies. However, our results show that quality of 
the movies is more important than quantity. Moreover, the exhibitor should prefer the 



high quality movies that are released early in the horizon so that they could schedule 
them more profitably during the planning horizon. This in turn suggests that the 
distributors of the strong movies should aim at releasing their movies earlier in the 
season so that the exhibitor has an incentive to play them longer. 


3.4 Comparison with Heuristics 

In this section, we compare the MDP approach with two benchmark heuristics. The 
first heuristic is motivated by theoretical results on structured policies (Puterman 
1994, p. 103) and the second heuristic is based on the relationship dilemmas in 
distribution channels (Swami, Eliashberg, and Weinberg 1999). 

3.4.1 Rank-based Optimal Policies and Greedy Heuristics 
Establishing the existence of optimal policies with a special structure is one of the 
principal uses of MDP methods. Such structured policies are useful because they are 
easy to implement, facilitate computation and allow investigation of the effect of 
model parameters. Examples of policies with simple structure include {s,S) policies in 
inventory control models and control limit or critical number policies in queuing 
control or equipment replacement models. Because of its similarity with the machine 
replacement problem, we expect that that an optimal policy for the movie replacement 
problem would have the form: In a given week, continue currently playing a movie if 
either it is in its ohligation period or its rank is below ' ' a critical number, otherwise 
replace it by a new movie. While we do not establish the optimality of this policy 
theoretically, numerically derived optimal policies suggest that this might be close to 
optimal. 

We investigate the performance of this policy in the following representative 
problem of eight weeks, five ranks and two movies (of the same type). The first 
movie, a (Type I movie), is playing at the beginning of the horizon, and the other one, 
h (Type 11 movie), is released at Week 3. The optimal movie replacement policy is 
presented graphically in Table 3.5. 


" Recall dial a higher rank miniber ofa movie indicates a wor.so po.silion. 



Table 3.5: Results of Rank-Based Optimal Policy Analysis 


Movie currently playing = a ('Fype 1) 

Movie available for replacement = b (Type 2) 

Actions 

U Retain a 
Re|)lace a 


Optimal Actions Space 


Decision epoch 
Number of weeks 
Movie a played* 
Ranks of Movie 

-,i,r 



'/ 2 

5 


3 ’ 4 7 5- ■ 

Oi 1 2 

1 2 3 4‘j' 'l T 4 5 12 3 Ts 





Decision epoch 6 ' 7 8 

Number of weeks 

Movie aplayed* 3 4 5 

Ranks of movie 

a\b 1 2 3 4 5 1 2 3 4 5 1 2 3 4 5 

V _ 

2 ' ' ■ ^ ' ■ 



* “ Beyond its obligation period of two weeks. 

The figure shows that in decision epochs 3-5 if Movie a is in Rank 1 or 2, then it is 
optimal to continue playing it no matter which rank Movie b opens in and in decision 
epochs 6-8 it is optimal to continue playing movie a in ranks 1-3. At decision epoch 3 
the optimal policy replaces a if its rank is lower than that of b. At subsequent 
decision epochs, the optimal policy departs slightly from this form. 

Similar results for similar other test problems motivate using the following 


variant of a critical rank heuristic: 



Rank Based Heuristic. Each week, if a new movie becomes available, replace the 
existing movie only if it is not in its obligation period and either its rank is above a 
critical value or the rank of the new movie is better than that of the one being 
currently playing.'^ 

Swami, Eliashberg, and Weinberg (1999) suggest that exhibitors sometimes 
face a channel relationship dilemma. Distributors want to find screen space for their 
new movies and pressure the exhibitor to book the latest release. Such an exhibitor 
uses a decision rule such as: 

Distributor Pressure Heuristic: Each week, if a new movie becomes available, to 
maintain a good relationship with the distributor. I will replace my existing movie, if 
it is not in its obligation period, with the new movie. If more than one new movie 
become available in a week, then I will replace the current movie by the highest 
ranked movie. Notice that the manager following this heuristics tries to accommodate 
the new movies released by all the distributors as long as he/she has “space” to show 
them. However, with some element of smartness in the decisions, the exhibitor 
replaces only by the best movies (in terms of ranks) available at any decision point. 

3.4.2 Analysis Results 

The comparison of MDP approach with rank based and distributor pressure heuristics 
was performed on twenty-five 16-movie problems. The expected values from using 
the heuristics was calculated using backwards induction with a fixed decision rule 
specified by the heuristic. Figure 3.3 shows the cumulative expected values of the 
MDP model and the two heuristics. 


Swami (1998) discusses a version of this heuristic and compares its performance to other decision- 
making approaches including a benchmark optimal policy. 



Figure 3.3: Comparison of MDP Optimal Policy with Heuristic Policies 



Problem Number 


‘ ■ MDP — ■ — Rank-Based Heuristic ■ A - Distributors' Pressure Heuristic 

Expected total revenues for three different policies over 25 problems. Mean for MDP optimal policy is 
2170, mean for Rank-Based Heuristic is 1842 and mean for Distributors’ Pressure Heuristic is 1751. 
Observe that the MDP optimal policy dominates the other two in all problem eases. 

As can be expected, the mean value of the rank-based heuristic (1842) is 
higher than that of the distributors’ pressure heuristic (1751). This is because in the 
rank-based heuristic, the exhibitor does not “accommodate” the distributors. This loss 
in value can be characterized as a “cost” of relationship management. However, the 
rank-based heuristic performs considerably worse than the optimal policy (2170) for 
the moderate size problems considered in this analysis. These results points towards 
the significant advantage that can be gained by using the MDP approach in larger, 
ongoing decision-making situations because the MDP approach captures more of the 
subtleties of a dynamic decision situation than simple heuristics. 


3.5 Implementation of the MDP Approach 

In this section, we discuss how the MDP approach may be applied in a realistic 
setting. To illustrate this application, we use the data of Swami, Eliashberg, and 
Weinberg (1999), which was publicly available in Variety (1989). The theater 
considered is the 84"’ St. Sixplex, a six-screen theater in New York. 

The movie scheduling recommendations are developed for the 1989 eight- 
week summer season. Ihe data available are ex post box-office revenue data, which 
had to be recast for the MDP approach in the current application. For this purpose, we 
considered 27 weeks of the year 1989 and a corresponding period from its preceding 
year 1988. We classified the movies of the two years separately in movie types 
following Jedidi, Krider, and Weinberg (1998). First, data were median split on the 
basis of the opening strength of movies. Then, each of these groups were median split 
on the basis of their decay rates. Thus, we obtained four “types” of movies on the 
basis of their opening and decay rates for both the years 1988 and 1989. Their 
corresponding expected revenues, initial and transition probabilities were generated in 
a similar manner as that followed in the simulation discussed in Section 3.3.1. 

To summarize, the transition matrices for the four clusters were estimated 
using 1988 data. Then, prior to the 1989 season, it is assumed that the exhibitor would 
know the list of the available movies, their release dates, the cluster that each movie 
fits in, and the corresponding movie parameters (e.g., probability matrices). However, 
the exhibitor does not know a priori which rank each movie would actually achieve in 
a given week.''* 

Applying the MDP model to this problem, we found that in six out of the eight 
weeks, the model selected the movie that had the highest revenue for that week as 
recommended by the integer linear programming based model SilverScreener of 
Swami, Eliashberg, and Weinberg (1999); in the other two weeks the MDP model 
selected the movie that had the second highest revenue. These results, though only 
illustrative and not conclusive, are encouraging because Swami, Eliashberg and 
Weinberg’s (1999) results partially use ex post box-office revenue data. The current 

This approach is based on Swami, Eliashberg and Weinberg (1999) and is also followed in practice 
by several exhibitors. 



results indicate that the MDP approach could mimic the “best” decisions in movie 
scheduling, where best is delined in terms of optimal nu>vie scheduling done using cx 
post data. Also, we presented the application for a single eight-week time window. 
Extensions of the approach for multiple screens and in a “adaptive scheduling” 
manner (Swami, Eliashberg and Weinberg 1999) could provide solutions in more 
realistic settings. 

3.6 Conclusion, Limitations and Future Research 

1 liis chapter addresses the need to do more integrative research between marketing 
and operations research and applies tools from both fields to solve an interesting 
management problem. We propose a modeling approach to analyze a retail 
management problem regarding the replacement of motion pictures on theater 
screens. The problem is complex because it involves stochastic elements in a dynamic 
decision making setting. One of the strengths of the MDP formulation is that it helps 
represent this complex problem by breaking it into separate components, which are 
easier to model. The resulting model provides theoretical and conceptual insights into 
the replacement phenomenon of perishable products such as movies. It then addresses 
some typical decision situations faced by movie exhibitors and provides optimal 
strategies for them to follow. 

The quality, quantity and availability analysis suggests that when there are 
more high quality films, the exhibitor’s screen (shelf-space) becomes more valuable. 
On the other hand, availability of additional low quality movies is not as valuable. 

We also investigated the possibility of characterizing the structure of rank- 
based optimal policies as a control limit policy. Our policy results depend on the state 
of both the current and replacement movies. This creates a two (or multi-)- 
dimensional decision space whose area depends on the relative attractiveness of the 
current and replacement movies, and the stage of the planning horizon. I'hese results 
suggest an easily implementable heuristic. But our evaluation of it shows that the 
optimal policy from the MDP yields significantly higher revenue. 

We now discuss some limitations of our study. We have assumed that rank 
transition probabilities are known and stationary. It would be interesting in future 



research to examine how learning and dynamic updating of probabilities would affect 
optimal policies. 


Another issue concerns the end-of-horizon effects. The end-of-horizon effects 
are usually summarized by an appropriate salvage value. Though we use zero salvage 
value in our example problems, the methodology allows any reasonable terminal 


value, h'or example we 


might replace zero by the expected total revenue of the movie 


playing in the last period. 


Finally, we restricted attention to a single screen case. However, interesting 
issues arise in case of a multiple-screen Ihcaler. Sonic oflhcse issues are allocation of 
different movies on different capacity screens, switching of the movies between 
screens, multiple screenings of the same movie, and so on. These are interesting 
research questions, which we have partially addressed in the next chapter. In addition, 
other industries pose similar imeresting problems involving the allocation of selling 
spQct to pciishiiblc pioducts with liniitcd lifetimes. 



CHAPTER 4 

MULTIPLE SCREEN PROBLEM: 
GENETIC ALGORITHM APPROACH 

4.1 Introduction 

I'he previous chapter dealt with a single screen problem. In this chapter we consider a 
ditlercnt capacity multiple-screen problem. Most of the movie theatres in the major 
cities in the world are multiplexes, that is, they have more than one screening rooms. 
The exhibitor is interested not only in the optimal use of any one of his/her screens 
but has to optimize the use of all the screens. We propose to address these issues in 
this chapter. We present use of Genetic Algorithm based heuristic in solving the 
multiple-screen problem, for which simple mathematical programming approaches do 
not provide easily tractable exact solutions. Some interesting issues arise in case of a 
multiple-.screcn theater. Since the capacity for each screen is different, the decision of 
whether or not to schedule a movie, is accompanied by a decision of which movie 
should be allocated to which screen. If only demand of a movie is considered, the 
allotment may not be optimal because of some other factors like contract term 
between the exhibitor and the distributor. Again, with the varying demand for each 
movie and a varying percentage of total revenue that the exhibitor will retain, the 
movies may have to be switched between screens. The issues like multiple screening 
of movies (i.e., the same movie displayed on more than one screen in a given week) 
are other interesting issues arising in this case. 

The chapter is organized as follows. In the next section we discuss the 
previous work specific to the problem. In Section 4.3, we explain in detail the exact 
problem and brielly explain the SilvcrScrecncr model developed by Swami, 
Eliashbcrg, and Weinberg (1999) which deals with a similar problem. We also explain 
the limitation of this model with the multiple-screen formulation proposed by Saxena 
(2000). We explain by an illustrative example, the situations in which the heuristic 
proposed by Saxena (2000) might give inferior solutions and justify the use of some 
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other heuristic. Towards this end, in Section 4.5, we explain the proposed Genetic 
Algorithm (GA) based heuristic. In Section 4.6, we discuss the computational study 
and the results obtained. In Section 4.7, we discuss the managerial implications and 
conclusions of our study. We also discuss the directions for future research. 

4.2 Literature Review 

In this section we discuss the literature that is very specific to this problem. Swami, 
Eliashberg and Weinberg (1999) deal with the problem of equal capacity multiple- 
screen theater scheduling. They propose an integer programming approach to this 
problem. Saxena (2000) extends their work by implementing the SilverScreener at a 
theater chain in The Netherlands. Saxena (2000) also addresses the limitation of 
model, which assumes equal capacity of all screens and proposes a conceptual non- 
linear model to relax the assumption. To address the issue of unequal screen 
capacities, Saxena (2000) proposes a screen allotment heuristic. We propose a GA 
based heuristic for this problem and compare the results of our heuristic with 
Saxena’s (2000) screen allotment heuristic. 

4.3 Problem Formulation 

4.3.1 The Exhibitor Problem 

In this section, we present the space-allocation problem faced by exhibitors, as 
proposed by Swami, Eliashberg, and Weinberg (1999). 'I'he above researchers 
presented the exhibitor's problem and its formulation in the North American context. 
Swami, Eliashberg, and Weinberg (1999) present the exhibitor's problem as follows. 
Every week, motion picture exhibitors have to make an important decision regarding 
the replacement of the movies playing at the screens in their theatres. The dynamic 
environment thus induced gives rise to the notions of decay and aging of movies. 
Decay is the intrinsic weekly decline in the box-office attraction and gross revenues 
(grosses in industry jargon) of a movie playing at a theater (Kridcr and Weinberg 
1998). Aging is the decline in the value, that is, gross generating power, of a movie 
from an exhibitor’s perspective if there is a delay (by week) in exhibiting the movie at 
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the theater. Aging, therefore, results in an opportunity cost of not being able to play a 
particular movie. 

1 he above scenario is further complicated by the nature of the contract 
between the distributor and the exhibitor. The basic structure of the contract is fairly 
standaid between dillerent distributor-exhibitor pairs although the individual terms 
may vary depending on the relationship between the two parties. A typical exhibition 
contract states a fixed obligation period and a differential revenue sharing scheme in 
different weeks between the distributor and the exhibitor. The obligation period 
limits the ability of an exhibitor to replace a movie with less than satisfactory 
box-office perfonnance in the initial weeks after its release. 

In a given week, the revenue sharing scheme splits the gross of a movie 
between a distributor and exhibitor by one olThe two rules: a) 90%/10% over house 
nut' with minimum gross percentage, or b) a movie type based contract term where 
the percentage revenue share for distributors and exhibitors depend on the type of 
movie (Swami, Lee and Weinberg, 2000). If the 90%/70% over house nut rule 
operates, then the distributor receives 90% of the gross after the exhibitor has 
deducted and retained the house nut amount. Accordingly, under this rule the 
exhibitor keeps 10% of the gross over house nut plus the house nut amount. The 
exhibition contract also contains minimum percentage figures as specified by the 
distributor for every week of the expected play length of a movie. These figures will 
be used if the minimum gross percentage sharing rule is invoked for revenue sharing. 
Under this rule, the whole gross box-oflice revenue amount (without house nut 
deduction) is split according to the specified percentage for that week. The splitting 
temis (in favor of distributor and exhibitor, respectively) specified by the distributor 
under a typical contract may appear as follows (see, for example. Squire (1992), p. 
315, for the movie Batman). 

90%/ 10% over approved house allowance with minimums of: 

First week at lOVol^OVo 
Next two weeks at 60%/40% 

' House nut is a small, negotiated amount, which the exhibitor receives from the distributor. It 
does not necessarily bear any relationship to the theater's actual expenses, and is only meant to allow 
for some cu.shion in the exhibitor's profit margins. 
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Next week at 50%/50% 

Next week at 40%/6()% 

Balance at 35%/65% 

The manner in which the box-office grosses are split favors the studio 
(distributors) in the first few weeks of the movie playing, but shifts to the exhibitor’s 
tavor later on. Distributors thus have a strong incentive to promote the movies 
intensively in their initial play period. On the other hand, the longer the exhibitor 
plays the movie, the larger becomes his/her share of the box-office receipts. At the 
same time, theater attendance for a movie typically declines the longer the movie 
plays. 

In addition to the revenue earned from the box-office gross, the exhibitor also 
earns some income from concession sales such as popcorn, candies, and soft drinks. 
The concession sales, however, depend on individual demands generated by the 
movies playing at the theater. The exhibitor does not share the concession income 
with the distributor. Most theaters have multiple movies playing at their multiple 
screens. Given the complexity of the revenue sharing scheme and the dynamics of 
movie availability and decision making, the managers ol'such multiplexes are faced 
with non-trivial decisions of selecting and scheduling movies on different screens in a 
fixed planning horizon. In the next sub-section, we discuss a mathematical model 
SilverScreener proposed by Swami, Eliashberg, and Weinberg (1999). 

4.3.2 SilverScreener: The Screens Management Model 
4.3.2. 1 Assumptions 

Swami, Eliashberg, and Weinberg (1999), in order to simplify the exposition, make 
the following assumptions. 

1. The availability of the movies to be released during the planning horizon is 
known in advance, 

2. 'I'he weekly revenues to be generated by the movies considered during the 
planning horizon can be estimated in advance, 

3. The replacement decisions are made on a weekly basis, 

4. All the screens in the multiplex are of equal capacity. 
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5. 


I here is no time lag between placing an order for a new movie and its arrival, 


Assumptions 1, 3 and 5 are not limiting and, in fact reflect current industry practice. 
Assumption 2 is a strong assumption about a priori knowledge of movie revenues. 
However, data collected from Variety suggest that managers have a reasonable 
estimate of box-office gross revenue of a movie. Equal capacity screens are assumed 
(Assumption 4), which addresses the questions regarding which movies to pick and 
how long to show them. In the model implementation section, we discuss Saxena’s 
(2000) approach to this assumption. We also explain the serious limitation in 
Saxena's (2000) approach and propose a GA-based heuristic to overcome it. 

4.3.2. 2 Definition of Variables 
W - length of planning horizon, 

H - number of screens in the multiplex, 

N - total number of movies considered during the planning horizon, 

Kj - Release date of movie j, 

dj - due date (if applicable) of movie y, 

- binary (decision) 0-1 variable which takes value 1 if movie j is scheduled for 

i weeks beyond its obligation period starting in week w, 

Py,„, - profit received by the exhibitor if is equal to 1, 

GROSS - box-office gross revenue generated by movie j in week w, 

POPjw - concession profit generated by movie / in week u’, 

VCjw - variable cost due to movie j in week w, 

FCw - fixed cost of multiplex in week vv, 

EXSHAREjw - exhibitor’s share of box-office revenue for movie / in week rv, 

OPDj - obligation period of moviey. 






C 


- house nut. 



4.3.2.3 The Model 
Tiino-liuleved Forimitatiuii 

We now discuss time-indexed formulation (Sousa and Wolsey 1992; Williams 1997) 
that is particularly useful for solving the exhibitor problem. This formulation is based 

on the idea of dividing the planning horizon [0 W] into W discrete intervals of unit 

length. To model the screen management problem, a binary variable is defined, 

which equals 1 if movie _/ is shown for i weeks beyond its obligation period starting in 
week vv of the planning horizon, and 0 otherwise. Notice that the obligation period 
constraint is included in the definition ol itself, ('or example, if obligation period 

of Movie Number 3 is 2 weeks, then .V301 = 1 implies that it is shown for 2 weeks 
starting in week 1 . 

The time-indexed fomiulation highlights some key differences between the 
exhibitor problem and typical machine scheduling problems. First, all movies do not 
have to be played in the exhibitor problem, while all jobs have to be scheduled in the 
machine scheduling problems. Accordingly, the proposed formulation is aimed at 
helping the exhibitor decide about two critical decision variables: choice of which 
movies and deciding play lengths of the chosen movies. Second, the following types 
of decision-making goals are found to be prevalent in machine scheduling problems: 
turnaround, timeliness, and throughput. In contrast, the exhibitor problem offers a 
situation in which the vseheduling decisions directly affect profitability, a more 
relevant decision making goal, both by affecting gross revenue and concession profits. 
Finally, the main parameters of interest in machine scheduling are the lengths of the 
Jobs, and their release and due dates. I'he exhibitor problem, on the other hand, deals 
with additional variables such as complicated contract terms and an exponentially 
decaying demand function. 
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Profit Function 

Pjiyy ^ total protit the exhibitor receives corresponding to each , is defined as: 

C = Z(- + POP^. - PC,J + /„_ * {O.l ♦ (GROSS,, -C) + n + 

(I exshare. „ * GROSS.,, , (4.3 . 1 ) 

j = i = 0,...,kj, w = r.,-,dj-SCRj.+\. 

where Q is a logical condition given by 

0;„, = (0.9 ♦ {GROSS p, - C) > (1 - EXSHARE j^;) * GROSS jJ, 
and 

lx = 1 i(X=TRUE 
= 0, otherwise, 

^7 “ * 0 ■ OPDj +1 = maximum possible number of weeks movie j can be 

shown beyond its obligation period starting in r, or any feasible week thereafter, 

SCRji = OPDj + i = total screening period for movie j if it is shown for i weeks 
beyond its obligation period, where i = 0, kj 

The profit expression proposed by Swami, Eliashberg, and Weinberg (1999) 
consists of two portions corresponding to the two different conditions of the contract 
temis^. The two conditions, operationalized by the logical variable /, follow directly 
from the revenue sharing tenns described earlier. In addition, there are variables for 
fixed and variable costs. Fixed costs (e.g., rent, lights, and weekly maintenance) are 
the costs that the exhibitor has to incur irrespective of the traffic generated by the 
movies playing in a particular week. Variable costs, such as salaries of the temporary 
staff hired for a particular movie, vary by movies as well as weeks. 

The above operationalization of Pjiw^' simplifies the solution procedure 

considerably since they can now be computed independently of the optimization 
routine. Further, the variables kj and SCRjj help us “cover” all feasible (P,x) pairs 
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for a movie. For example, suppose movie ; has parameters: OPD^ =2, r =1, and 
d j - 4. It the movie is scheduled in Week 1 (i.e., w =1) , then it can be shown for 2, 3 

or 4 weeks, that is, for 0, 1 or 2 weeks beyond the obligation period. Since kj =2, i 
varies from 0 to 2. The corresponding values of 60^,- (= 2 + /) are 2, 3, and 4, 

respectively. Therefore, the corresponding {P.x) pairs are {Pjoi, xjoi), {Pjn, xju), and 
{P/ 21 . respectively. 


Demand Function 

Swami, Eliashberg and Weinberg (1999), Krider and Weinberg (1998), and Sawhney 
and Eliashberg (1996) use an exponentially declining demand curve to estimate the 
number of visitors attracted by a movie. 

Demand (number of visitors) = F7S'/rO/?/„= (4.3.2) 

where a,- > 0 and jij < 0 are opening and decay factors for Movie y, 

and e ~ normal (0, a^). 


Problem Statement 

The time-indexed formulation of the exhibitor problem is presented as follows: 

maxY Z Z PjiwXjhv (4.3.3) 

/ = ! /=0 w-f'i 

kj (lj-SCRji+\ 

subject to Z Z = (4.3.4) 

/=0 »■=;■; 

SS t (4.3.5) 

/ I i 0 ir ll 


< (jj < d. - SCR„ + 1, j = 1,. ..,N; i = 0,...,kj (4.3.6) 

^ The SilverScreencr does not consider the movie based contract term, however while comparing the 
performance of proposed GA-based heuristic with Saxena’s (2000) heuristic, which uses 
SilverScreencr in its first step, we made suitable modifications in the SilverScreencr model to consider 
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■^‘/nr {0^1} (4.3 7) 

In the abo^e model. Statement (4.3.3) denotes the objective function, which is to 
maximize cumulative profit o\ er the season. Constraint (4.3.4) ensures that a movie is 
plav ed in only consecutive weeks. It also allows a movie not to be scheduled at all. 
The next constraint restricts the total number of movies, scheduled in any week of the 
planning horizon, to the total number of screens in the multiplex. In doing so, it sums 
up all the movies, which are released earlier than or in the week under consideration. 
The set of inequalities denoted by Equation (4.3.6) is an indexing constraint, which 
restricts the variable q, in Constraint (4.3.5) to feasible values. Constraint (4.3,7) 
defines be a binary variable. 

4 J.2.4 SilverScreener Implementation at Pathe, The Netherlands 
Screen Capacity 

In their formulation Swami. Eliashberg and Weinberg (1999) assumed all the screens 
of equal capacity. However in the implementation of the model. Sa.xena (2()()()) 
follow ed the follow ing simple heuristic for allocating the movies on different capacity 
screens. While evaluating the performance of GA based heuristic, we will compare its 
performance with Sa.xena's (2000) heuristic. 

Screen Allotment Heuristic (Saxena 2000) 

The statement of the heuristic is; 

In a particular week, allocate the movie with highest number of visitors to the highest 
capacity screen, the movie with the next highest number of visitors to the next highest 
capacity screen, and so on. 

In other w ords, the application of the model first chooses a set of movies, and then the 
mo\ ies are allocated to tlie tlieater screens in tlie order of tlteir capacities, live screen 
allocation heuristic is similar to tlie managerial decision making in some instances. 


the mo\’ie Ixiscd contract tenns. 



4.3.2.5 Illustrative Example to Show Limitation of the Heuristic by Saxena (2000) 

We illustrate the limitation of the screen allotment heuristic proposed by Saxena 
(2000) with a simple example. Consider a two week two screen problem. Suppose 
that there are three mo\ies a\ailable lor playing Ihc demand for these mo\ies for 
two weeks is as shown below. 


Demand for Movie 

Week\Movie 1 2 

1 ~ ITtK) ' ' " ^ 1200 

2 200 ' 1000 


3 

500 

300 


The exhibitor share in net revenue for each movie in each week is assumed as shown 
below. 

Percentage Exhibitor Share 

”1 ' "l ' 3 

10 ’ 15 ’ ~'30 

30 20 ' 50 


Runlcngth\Mo\ ie 

' T 

2 


Let the screen capacities of the tw o screens be 700 and 200, respectively. 

The SilverScreener will select Movies 1 and 2 for the first week and Movies 2 and 3 
for the second week. The Screen allotment heuristic proposed by Saxena (2000) will 
allocate Screen 1 to Movie 1 in first week and to Movie 2 in second week. Similarly 
Screen 2 will be allocated to Movie 2 in Week 1 and Movie 3 in Week 2. The revenue 
earned will be as shown below. 
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Slot 

Movie 

* Revenue 

Week 1 Screen 1 

1 

700*0.1=70 

Week i - Screen 2 

2 

200*6.15=30 

Week 2 Screen 1 

2 

700*0.2=140 

Week 2 - Screen 2 

3 

200*0.3=60 

Thus the total revenue for exhibitor will be 300. 
follows 

A better schedule can be generated as 

Slot 

Movie 

Revenue 

Week 1 - Screen 1 

3 

500*0.3=150 

Week 1 - Screen 2 

2 

200*0.15=30 

Week 2 - Screen 1 

2 

700*0.2= i 40 

Week 2 - Screen 2 

3 

200*0.5=100 

This schedule will generate the exhibitor revenue of 420. Thus, the example illustrates 

the limitation of screen allotment heuristic as proposed by Saxena (2000). Since GA- 

based heuristics are search procedures we 

might reach 

the better solution by 


employing such techniques. 

4.4 Multiple Screen Formulation (Saxena 2000) 

SilverScreener model assumes all the screens of equal capacity (Assumption 4) for 
simplicity. The screen capacity is, of course, a critical variable to the problem 
formulation only if the demand of a movie exceeds screen capacity. Indeed, the 
primary reason for the exhibitor to have dilferent screen sizes is that the attendance 
varies for different movies and having different sizes increases scheduling flexibility 
(while saving capital cost in the original construction of the theater), A movie may be 
switched from a higher to a lower capacity screen'. We now discuss the extension to 
the model as proposed by Saxeiia (2000) for an unequal capacity case. 
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Derinition of Variables 

T - l.cngtli of planning horizon, 

H - Number of screens in the multiplex, 

hJ - Number ot movies considered during the planning horizon, 

r, - Release date of movie j 

d j - Due date of movie j 

cap^ - Capacity of screen s 

- Binary (decision) 0-1 variable which takes value 1 if movie./ in scheduled 
on screen .v in week u 

y ■ Binary (decision) 0-1 variable which takes value 1 if movie / in scheduled 

in week w 

- Profit received by the exhibitor if jc is equal to 1 
T’/,v//o/-[j,u] - number of visitors for movie j in week u 
JJoor[jAi] - distributor's share in the box-office revenues for movie j in week u 

Problem Formulation 

Saxena (2000) defines , the exhibitor receives corresponding to , as 
Pjsu = [Min { ’ visitor [j, u]}] * (2 + 13.5*0.89725*(l-floor [j, u])) 

Since the implementation of SilverScreener by Saxena (2000) was in context of The 
Netherlands, the profit function is according to the practices of The Netherlands’ 
motion picture industry. The above problem is formulated as an integer-programming 
problem. I'hc tinic-indc,\od formulation is presented below. 

NUT 

Objective function: Max 

/=:1 55=1 uM 

Subject to: 

1 V J.u (4.4.2) 

5 = 1 


^ The opposite case may also occur when a movie’s demand builds on word-of-mouth in the later 
weeks though rarely in practice. 
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.V // 

II- 

I“1 S”1 

V 

u = I..T 

(4.4.3) 

Jsu ~ 

0 V 

j.sand d.<u<rj 

(4.4.4) 

H 

s~\ 

1 < Max 1 1 

^ " 1 
“Zt/u. >Z^./- f ^ 

U-! J s=\ J 

-T-l.j (4.4.5) 

th. 

1/ II 

- II-, 

V=1 5=1 

V j,u=L.T 

(4.4.6) 

N 

= 1 

V s, u 

(4.4.7) 


> e(0,l) 


(4.4.8) 


In the above model, Constraint (4.4.2) ensures that in a particular week u any movie j 
can only be scheduled at most on one screen, that is, it cannot be scheduled on more 
then one screen in a given week. The inequality sign is to ensure that the movie may 
or may not be scheduled. For example if x,,, =1 i.e. movie 1 is scheduled on Screen 1 
in Week 1 then the Constraint (4.4.2) +x, 2 , +.r, 3 , +x^^^+x^^^ +a:i 6, <1 implies 

^121 = ■’^ 1.21 “ -^141 = ^151 = -*161 = 0 as these are binary variables and can take the values 
Oor 1. 

Constraint (4.4.3) restricts the total number of movies scheduled in a given 
week to the total number of screens in the theater. Constraint (4.4.4) simply ensures 
that the movie can not be scheduled before its release date and after its due date. In 
most instances the movie runs in a continuous manner and it is very rare that the 
movie mns for few weeks and then resumes after a break of one or more then one 
weeks. 

Constraint (4.4.5) deals with the continuity of the movie mn. It restricts the 
movie to play in breaks. So this constraint schedule the movie in such a manner that 
the movie run is continuous. I’here are two situations that may arise. We discuss these 
situations by taking a simple example: 
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Case 1 ; Movie plays for few weeks and then a break of one or more w eeks 
Suppose the release date of' the Movie j is Week 1 and it runs in Week 1 and 2 and 
does not run in Week 3. 


So = h.v,, = = 1 and y’,,=0 

So due to constraint (4.4.5), 


■V,-, 4 +-^ 24 +-^',;, 

34 +'^744 +-^'754 +-^764 ^ 

as movie is not played in Week 3, therefore, 

''-■7.3 + •'■ 7'23 + -^73 

3 •^*y43 •^./53 -^763 ”” ^ 

Hence, 



34 +'>'’,/44 +'^';54 + '^;;64 ^ { 

•^'714 + ■'fy :4 + -’^7 

34 '^ /44 *^,/54 ^ /64 ~ 0 


(l-l-l-O), 1 

JU + ^;43 ■*" -^753 '^y63 


Hence the movie j cannot he scheduled in Week 4. 

Case 2: Movie not scheduled on its release date 
Suppose movie is not scheduled in Week 1, 2 and 3. 

0 = ^ >';l = yj2 = >’y3 = ^ 

therefore, according to constraint (4.4.5), 

(l-O-O-O), 


'^■.,14 + Xj2A + Xj24 + '"^244 + ■’Cy54 + ^>64 ^ 


•^713 ■^/23 -’'•,/33 ^/43 -’^,;53 •^>63 


As movie is not played in Week 3, therefore 


•’^713 

+ •’^723 + ^233 + ^743 

+ ^753 

+ - 0 

So 

•' ;14 4 •' ,24 4 •' , 34 

^ -^44 

f .V,54 + .V,„4 


< 1 

Hence the Movie j can be scheduled in Week 4. 

However, Constraint (4.4.5), introduces nonlinearity into the system and 
hence, difficult to solve exactly, fhis justifies the use of a heuristic method for 
addressing this problem. 
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4.5 Genetic Algorithms 

The concept of Genetic Algorithms (GA) was first proposed by Holland ( 1975 ). The 
basis for the algorithm was the observation that a combination of reproduction and 
natural selection allows nature to develop living species that are highly adapted to 
their environment. To apply GA to a problem, the candidate solutions have to be 
formulated in teirns of strings. 

Genetic Algorithm is basically a search algorithm in which we search for an 
optimal solution point for a certain objective. At each iteration, known as a 
generation, we have a fixed number of solution points. Each solution point is called as 
a chromosome. The fixed number of chromosomes at any generation is called as the 
population size. With each generation wo move on to belter and better quality of 
chromosomes, with respect to the objective decided. The maximum number of 
generations is decided by certain convergence criterion. 

Initially we start with chromosomes that are generated either randomly or by 
certain heuristic. This group of initial chromosomes, which are equal in number to the 
population size, is called the initial population. Each chromosome is evaluated for the 
objective that has been decided. To move on to the next generation, we use certain 
genetic operators. We select some chromosomes from the current population, which 
are used for generation of the next generation. There are different ways to select these 
chromosomes. The basic objective of the selection operator is to increase the 
proportion of good chromosomes in the mating pool. A mating pool is a group of 
chromosomes, which will undergo certain operations to produce the next generation 
chromosomes which are better than the present ones, when evaluated against same 
criteria. 

The genetic operators used to produce the next generation which has ‘better’ 
chromosomes, are crossover and mutation. Each of the chromosomes of the mating 
pool does not undergo mating. For the crossover operation, each time a pair of 
chromosomes is selected and with a certain prefixed probability, called probability of 
crossover, they are allowed to crossover. In crossover, the selected chromosomes 
exchange elements with each other. The underlying principle of crossover is that with 
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the exchange of elements between chromosomes we may get better ones. A crossover 
site in a chromosome, along which the elements are exchanged, is selected randomly. 

1 he mutation operator is used for diversification. If only crossover operator is used 
the seatch gets tonimed to cettain space and solutions can converge at the local 
optima. In mutation we allow the chromosome to diversify from its parent. The 
mutation operator is generally applied bit-wise, that is, we have a pre-defined 
probability, called mutation probability, with which an element in the parent 
chromosome is allowed to mutate. The element takes a diverse value at random and 
thus results in the diversification of the population. The approach therefore operates in 
the following manner. 

4.5.1 Algorithm Process , 

The basic GA process can be visualized as follows (Michalewicz 1992, Grefenstette 
1986). Here POP{r) indicates population at generation t. 

Step 1 : t <— 0 

Step 2: Generate initial population, POP (t) either randomly or using a heuristic. 

Step 3; Evaluate each of the strings in POP(t ). 

Step 4: IF stopping condition satisfied THEN stop ELSE continue. 

Step 5 : / <— / + 1 . 

Step 6: Select strings from POP(/ -1) to create POP( t ). 

Step 7; Apply genetic operators on strings in POP (/ ). 

Step 8; Evaluate each of the strings in POP(t). 

Step 9: GO TO 4. 

I'he algorithm thus works on the fact that a “good” chromosome contains some 
features (sub-strings) that are desirable and hence contribute to its being evaluated 
high. With the exchange of genetic material between two good strings, the expectation 
is that the offspring generated will be good. Thus, as the iterative process continues, 
we move to better and better population of solutions. 

GAs have been applied in a variety of fields. Goldberg (1989) has used GA for 
ojitimization ol pipeline systems; Cleveland and Smith (1989) for scheduling flow 
shop releases; Liepins and Potter (1991) use GA for multiple fault diagnosis; 
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Balakrislman and Jacob (1996) use GA for product design; Deb, Chataborty and De 
(1999) use GA for model based object recognition; Joshi (2000) uses GA in search 
Space of hypercubes. 

4.5.2 GA in Multiple Screens Problem 

By the application of Genetic Algorithm to solve the multiple-screen problem, we 
develop an approach for finding a schedule to maximize the total revenue the 
exhibitor will obtain. 

For the multiple-screen problem, GA views sequences of movies to be 
displayed on screens for each week, for the total planning horizon, as a chromosome 
(the candidate solution), which in turn is a member of a population. Fitness of each 
chromosome is measured by the value of the objective function (i.e. the total revenue 
earned by the exhibitor.) The GA performance depends greatly on the choices of the 
different GA parameters, population size (pj, the probability of crossover {pc) and the 
probability of mutation (p,„). Selection of these parameters is itself an optimization 
problem and these parameterization have also been noted to be problem specific 
(Kumar, Bagchi and Sriskandarajah, 2000; Davis, 1991; Grefenstette, 1986;). We use 
a design of experiments (DOE) approach (Bagchi and Deb, 1996) to quickly identify 
the good settings for ps, pc, and p„, before using GA to solve our problem. 

We define an h x w vector to represent a chromosome for an n screen w week 
problem. For example, a chromosome for a two-screen two-week problem can be 
represented as nn ms mi m 3 ■ The chromosome represents that 
Movie mi is played in Weeks 1 and 2 on Screen 1, Movie ms is played in Week 1 on 
Screen 2 and Movie m 3 is played in Week 2 on Screen 2. The GA begins with a 
population (population size ps) of randomly generated vectors. The fitness of any 
chromosome is then evaluated based on the total revenue for that chromosome. To 
calculate the total revenue corresponding to a movie when placed in a particular 
week-screen slot, the number of visitors is taken as the minimum of the respective 
screen capacity and demand of the movie in that week. Using this number and the 
contract term, we calculate the revenue earned by that movie in that week. Total 
revenue is the cumulative of all such week-screen slot corresponding revenues. 
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screen slot corresponding revenues. The revenue that the exhibitor earns from a week 
1 / screen .v slot if movie / is scheduled in this slot is given by 
Re vemwj,,^ = [Mi,, {cap^. , demand }]*(!- floor . ^ ) 

wheic is pciccntagc shaie ol distributor lor movie /, when run for weeks r and 

Ciips is capacity of screen 5 and deinandju is the demand for movie _/ in week u. 

4.5.3 GA Parameters 

We now briefly provide overview of Genetic Algorithm parameters employed in the 
current study that are particularly germane to the multiple-screen problem. We begin 
by discussing the operationalization of GA operators to the specific class of 
chromosomes dcllned above. 

Selection: We use the method of tournament selection. In this type of selection, at a 
time two chromosomes are selected randomly from the current population and 
compared on the basis of their fitness value (total revenue). The one with high fitness 
value is added to the mating pool. This procedure is repeated for ps number of times, 
and each time the winner is added to the mating pool. Thus we maintain the total 
population size p, in the next generation. Since the selection is based on the fitness 
value the ‘good’ chromosomes have a higher probability of getting selected for the 
mating pool as compared to bad ones. This increases the proportion of good 
chromosomes in the mating pool. These selected members are given a chance to mate 
according to the probability of crossover (/i,-). 

Crossover: A key consideration in designing the crossover operator is avoidance of 
infeasible solutions. First, we select the mating pairs of chromosomes at random from 
the mating pool. Then with a probability pc we allow them for mating. Then, we use 
“one-point” crossover (Murata and Ishibuchi, 1994) in which one point is randomly 
selected for dividing one parent. This point is the crossover site. Due to the complex 
structure of the chromosome, we restrict the crossover site to be a starting point of any 
week. I'hus the crossover site can be either starting point of the first week or second 
week and so on. This avoids the in-feasibility of resulting chromosome due to release 
date constraint. The elements on the left-hand side of the crossover site are transferred 
from the parent to one child. The elements of right hand side for the child are the 



Figure 4.1(a) Parent Chromosomes before Crossover 
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With crossover point as starting of third week, the children formed will be as shown 
in Figure 4.1(b). 
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Week 


Figure 4.1(b) Children Chromosomes after Crossover 
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Figure 4.2 Mutation: Repiacement by a Movie Not Previousiy Scheduled. 
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Figure 4.3 Illustration of use of Repairing Function 
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4.5.4 Design of Kxpei iineiUs: Parameter Optimi/ation Approach 

The GA parameters, namely, the population size ps, the probability of crossover pc, 
and the probability of mutation p„„ have a great impact on the performance of the GA 
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based heuristic. Mere we explain the Desifru ofExpc 
three tiA paraineters. I'ach oi' these parameters 


rimenis approacli to optimize the 
has twt) settings Ibr which the 


perrormance of GA is tested. I he experimental setup is as shown in Table 4.1 


There are eight experiments that cover all combinations of settings for these 
three factors. I he ‘O' setting for a factor indicates that it is set to a lower amongst the 
two values, and the ‘ 1 ’ setting indicates that it is set to a higher value. For each 
experiment, the OA is run for a fixed number of generations and the revenue 
corresponding to the best schedule is noted down. Then for a factor, the sum of 
revenues corresponding to all the experiments in which the factor was set at ‘0’ is 
compared with that of experiments where it was set to ‘1’. The setting of the factor 
with higher value for the sum is then considered better, 'fhis is done for all the fiictors. 


Now if setting ‘1' is proved to be better for a lactor, we make it as a ‘0’ setting and 
make some higher value as a ‘ 1’ setting for the factor, thus moving in the direction of 
higher values. vSimilarly, if ‘0’ setting proves better, we move in the lower direction 
for that particular factor. With these new settings the process is repeated. If moving in 
any direction from a setting for a factor degrades the performance of GA, we 
conclude the setting for the factor as optimal. 


Tabic 4.1 Experimental Setup for GA Parameter Optimization 


Experiment 

Ps 

Pc 

Pm 

1 

0 

0 

0 

2 

1 0 

0 

1 

3 

0 

1 

0 

4 

■ ' 0 

1 

1 

5 

i 

0 

0 

... 
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0 


7 

’ 

i 

0 

8 


'i 

1 
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4.6 Perforiiiance Evaluation 

Saxcna (2000) reports that tlic integer programming model based heuristic worked 
quite well for mulliple-sercen problem. However, as illustrated in the illustrative 
example ot Section 4. 3. 2.4, his approach may lead to inferior solutions since first the 
movies to be scheduled are selected by the integer programming model, without 
considering the screen capacities and then the screens are allotted to the movies by 
considering only their demand. In allocating the screens to movies an important factor 
of contract temi between the exhibitor and distributor is neglected. 

'fhe GA approach that we propose is a search method, and hence may lead to actual 
optimal solution. I o dcnionstiate the eilectiveness of our approach, we need to show 
that under certain conditions of the input parameters, a GA-based heuiistic, which 
explicitly consiilcrs the capacity variations in the screens works belter. 

We now explain the metliodology to compare the performance of the proposed GA- 
based heuristic with the Saxena's (2()()()) heuristic, d'he problems examined for this 
purpose were generated in a similar manner of the previous chapter. Movies are 
classified into four “types" based on their opening strength and decay rate. As 
explained in the implementation section of the MDP model in previous chapter, the 
data of 84''' St. Sixplex, a six-screen theater in New York, was used to estimate the 
opening strength and the decay rate of each of these types of movies. We considered 
27 weeks of the year 1989 and a corresponding period from its preceding year 1988. 
We classified the movies of the two years separately in movie types following Jedidi, 
Krider, and Weinberg ( 1 998). First, data were median split on the basis of the opening 
strength of movies, 'riien, each of these groups was median split on the basis of their 
decay rates. The median values of the data for each of the four groups were used as 
the representative opening strength and decay rates for the groups. We estimated the 

demand pallcrn for each type of movie hy using llic e(]ualion where a 

represents the opening strength, represents the decay rate, and t represents the 
number of weeks after the release of movie. The constant multiplier 5 was chosen 
arbitrarily as we were interested in getting only the representative demand curve for 
each type of movie and not the actual demand for these movies types. 
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We first discuss the factors that we considered to be significant in deciding the 
iniproN cnicnt of (iA-bascd heuristic over Saxena’s (2000) heuristic. The constraint of 
capacity ol scieeus is the extra constraint that is considered by the GA-based 
heuristic. Therefore, we considered screen capacity to be one of the significant 
(actors. Second, we lound in our preliminary analysis that SilverScreener has a 
tendency to select Type I movies'* in their preliminary weeks because of their high 
opening strength. However with the added constraint of capacity, the value addition in 
the objective function corresponding to Type I movies might be reduced. Therefore, 
we select the number of Type I movies in a scenario as another factor to be analyzed. 

Ihe third lactor we analyze is the contract term between the exhibitor and the 


tlisli ibutoi , which decides the actual rew^enue earning strength ol the movies (or the 
exhibitoi. lo sunimatizc we study the impact of the variation in the capacity of 
screens, the numbei of 1 ype 1 movies released and the contract term between the 
exhibitoi and the distiibuloi, on the improvement shown by GA-based heuristic over 
Saxena’s (2()()()) heuristic. 

We examine two levels, high and low, of the capacity factor. To consider the 
typical distribution o( capacity over screens, we use the screen capacity data of a 
theater in The Netherlands. In the high capacity case, we assume eight shows in each 
screen room per week and in low capacity we assume four shows per week^. Table 
4.2 gives the capacities we assumed for each of the cases*’. 


Table 4.2 Capacity of Screens 
Capacity of Screen 


C'asc 

i 

2 

3 

4 

5 

6 

High Capacity 

3472 

2736 

1728 

1208 

1112 

904 

Low' Capacity 

1736 

1368 

864 

604 

556 

452 


'' I'ypc I movies arc the movies with high opening strength but decay rapidly. 

’ Again thc.se numbers of shows do not exactly represent the actual number of shows in a theater, and 
are scaled according to demand of movies in initial weeks. 

^ Note that in the two levels of Capacity we do not change the proportion of distribution of capacity 
over screens. Studying the impact of change in this proportion can be an interesting study and 
recommended for future research. 
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1 0 operationalize the proportion of 1 ype I movies in a scenario, we classify 
the sconuiio as a hi^h dci iiy [iropi’i ty scenario in which there are ten or more rapidly 
decaying movies (Type 1). Similarly, if there are three or less number of Type I 
movies, we classify the scenario as low ‘decay property’ scenmo. 

We consider the performance of GA under two types of contract terms that 
could be used, namely, first a 90%/10% contract term (common for all movies) and 
second movie type based contract terms. 

The 90%/ 10% contract term is used with the minimum percentage share as 
explained in the Section 4.3.1, with a house nut of 600 (scaled as per other factors. 
Squire, 1992). In their study of contract terms Swami, Lee and Weinberg (2000) 
suggest that movie type based contract terms result in improved channel coordination 
and hence recommend the same to the distributors. A typical movie based contract 
tenn considered by us is shown in Table 4.3. 


Table 4.3 Movie based Contract term 
Exhibitor’s % Share in Total Revenue for Movie Type 


Week 

I 

II 

III 

IV 

r 

10 

15 

25 

25 

i 

30 ‘^ ■ 

20 

40 

40 

3 onwards 

' '50 ~ 

~ 35 

50 

50 


These are according to the industry practice. Movies of Type I open with high 
demand and decay rapidly, so distributors of the.se movies have a tendency of 
retaining maximum share of revenue in their initial weeks. Since l ype 11 movies 
retain their strength for long, the distributor offers less share of revenue from these 
movies to the exhibitors. Movies of Type III and IV open with low strength and never 
have less demand, so to promote exhibitors to play these movies, distributors offer 
more share in the revenue for such movies. 

The input parameters required to generate problems for various simulation analyses 
are as follows: 
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i. length of planning horizon 

ii. number of screens at the theatre 

iii. quantity of movies (i.e., total number of movies considered) in a scenario 

iv. quality of movies 

V. release dates of movies 

To generate different simulation problems, we assign fixed values to the first 
three parameters, namely, the length of planning horizon, the number of screens in a 
theater and the quantity of movies released during the season, which are common to 
all problems. We set the length of the planning horizon to be 8 weeks (e.g., 
representing the two peak months of summer season) and the number of screens to be 
6. The quantity of movies released is also kept constant. We assume release of 4 
movies in each week. To begin with we assume that there are six movies already 
playing on the six screens at the end of the previous season and that the exhibitor can 
continue to play any of them in the new season.’ Therefore the movies released in 
week one of the season are not guaranteed of exhibition unless the exhibitor finds it 
more profitable to replace some of the currently playing movies with the newly 
released movies. The release dates of movies are thus fixed by assigning the number 
of movies released per week to 4. 

To randomize the assignment of a movie to a movie type, we use the relative 
proportions of Type I (19%), II (7%), III (38%), and IV (36%) movies in Jedidi, 
Krider, and Weinberg’s (1998) sample. Specifically, a random number is drawn 
among integers between 1 and 100 for each movie. Then the movie is assumed to be 
Type I if the random number is between 1 and 19, Type II if it’s between 20 and 26, 
Type 111 if it's between 27 and 64, and Type IV if it is between 65 and 100. 

4.6.1 Simulation Setup 

We study the effect of capacity of the screens (high/low), type of contract term 
(90%/ 10%/movie based), and decay property (high/low) for the period, on the 

’ 'riuis there are .t8 movies in caeli siimilatcd season. 


difference between exhibitor s cumulative revenue from proposed heuristic and 
Saxcna's (2000) heuristic. We use a complete 3 (treatment) X 2 (level) factorial 
design. Movie based contract term is denoted by type 1 and 90%/ 10% contract term is 
denoted by type 2. 

4.6.2 Methodology 

Using the above scheme of simulating problems, a total of 58 problems were 
generated. Based on the input data, the problems were assigned to their respective 
type of experiment. The GA model of Section 4.5 was coded in C++ language 
(Appendix 2) and run on an Intel Pentium III class computer. The value evaluations 
were done according to the cumulative revenue criterion presented in Section 4.5.2. 

4.6.3 Simulation Results 

Table 4.4 summarizes the simulation results in several ways. Table 4.4(a) gives the 
mean (and standard deviation) of improvement shown by GA-based heuristic over 
Saxcna’s (2000) heuristic, under each factor level combination and Table 4.4(b) gives 
the mean and standard deviation of improvement shown by GA-based heuristic over 
Saxena’s (2000) heuristic, for each main effect. The ANOVA results in Table 4.4(c) 
show that two interactions are not significant and that all three main effects (capacity, 
decay property and contract term) and an interaction of capacity and contract term are 
significant at the .05 level. 
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ruhlc 4.4 («): Mean liupioveinciU in ( uimiluUve Uevcuucs by (JA-based 
Heuristic over Saxena’s (2000) Heuristic for Experimental Setup 


Exp! 

Capacity^ 

Contract 

Decay 

Number 

Mean 

' % 



term^ 

Property' 

of 

improvement 

improv- 





problems 

in cumulative 

ement 






revenue* 


1 

0 

0 

0 

6 

4354(1183) 

42 

2 

1 

0 

0 

6 

3230(737) 

, 21 

3 

~o" 


0 ' 

6 

2488 (544) 

22 

^ 4 

0 

0 

1 

6 

5859 (449) 

56 

5 

1 

V 

0 

6 

2181(843) 

' 13 


0 

r~Y 

~ 1 

"6 

3695(551) 

35 

7 


0 

1 

6 

3595 (511) 

21 

8 


” 1 

1 

6 

3609 (397) 

23 


Total = 

48 

a - High (1) setting implies high capacity levels, otherwise low (0) setting 
b - Setting 0 implies movie based contract term, setting 1 implies 90/10 contract term 
c - High (1) setting implies more than 9 type I movies in the scenario, Low (0) setting implies less than 
4 type I movies. 

* Standard deviation values are shown in parentheses. 

♦♦ The total number of problems is 48 as 10 out of the randomly generated 58 problems 

could not be assigned exactly to an experiment type because the number of type I movies was in 

between 3 and 9. 



Table 4.4 (b): Mean Improvement in Cumulative Revenues by GA-based 
Heuristic over Saxena’s (2000) Heuristic for Different Settiii}»s of Parameters 


Factor 

Setting 

Mean Improvement in 



Cumulative Revenue* 

Capacity 

High 

3154 (622) 

' 

Low 

4099(681) 

Contract term 

90%/10% 

29 93 (583) 


Movie based 

4260 (720) 

Decay Property 

High 

4190 (477) 


Low 

3083 (826) 

* Standard deviation values arc shown 

in parentheses 



Table 4.4 (c): ANOVA Results for Improvement of GA-based Heuristic over 

Saxeiia’s (2000) Heuristic 


Source of 
Variation 

Sum of 
squares 

Degrees of 
freedom 

Mean square 

F() 

Capacity (A) 

10729534 

1 

10729534 

6.95‘ 

Contract term 

1 924320 i 

“l 

"19243201 

12.42’' 

m 





Decay 

'''15223521'" 

1 

15223521 

9.86 

Property 

if') 





AB 

”““6727518 

1 

6727518 

4.37 

AC 

632502 

i 

632502 ' 

0.4 

BC 

^438.536 

1 

"“137536 ' 

()l8 

ABC 

" 1386520 

”"””i 

1386520“ 

0.89 

Error 

18172536 

40 

45430 


Total 

'725^768 

■ 47 



• Significant at 5 percent. 




4.6.4 Discussion 




Intuition suggests that the factor of capacity should be the most important factor 


improvement of cumulative revenue by GA-based heuristic as compared to Saxena s 
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(2000) heuristic. The reason for this argument is that, the only difference in the 
Saxcna s (2000) heuristic and the proposed GA-based heuristic is that of a capacity 
constraint. However, our results show that other factors such as contract term and 
decay property of scenario are also important. In fact, the contract term is most 
significant factor followed by the factor of decay property and the capacity factor. The 
only interaction that is significant is the capacity and contract term interaction. In the 
90%/ 10% contract term, the exhibitor gets more or less an equal percentage of share 
for all types of movies. Therefore, in case of 90%/10% contract term, the only factor 
which remains important in optimal selection of movies is the demand for a movie. If 
the capacity constraint is applied in such situations, the Screen Allotment heuristic of 
Saxena (2000) might also provide good solutions. When the movie based contract 
term is applied, a movie of Type 1 which has a high demand at opening but, a low 
percentage of share for exhibitor in the initial weeks, does not remain optimal after 
the application of the capacity constraint. 'I'his is because the capacity constraint 
limits the demand for any type of movie to a maximum of screen capacity, while the 
percentage share of exhibitor in the total revenue, is high for other types of movies. 
Thus the fiictor of contract term between the exhibitor and distributor is most 
important factor. 

The second factor found important is the decay property of a season or the 
total number of rapidly decaying movies in a season. It is important to note that the 
setting for this lactor were very extreme settings. I'his may have introduced an 
artificial significance to this factor. Since Type I are the rapidly decaying movies the 
distributor of such movies have the tendency of retaining maximum revenue, from the 
first few weeks exhibition of these movies. So the initial exhibitor share for such 
movies is less. Each movie of Type I scheduled by Saxena’s (2000) heuristic leads to 
degradation of the schedule, because first, the exhibitor share in profit is less and 
second, these movies are allotted the screens with maximum capacities because of 
their high demand. Therefore, the number ot such movies in a season affects the 
improvement shown by GA-based heuristic over Saxena’s (2000) heuristic. 

The only interaction that was significant was the capacity and contract term 
interaction. The proposed GA-based heuristic shows more improvement over 
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Saxena’s (2000) heuristic if this interaction is high. When the capacity is set at low 
value (0) and the contiact tcini lactor is also set at movie based contract term (0), then 
their interaction is high. This combination of settings has a very large improvement 
shown by GA-based heuristic over Saxena’s (2000) heuristic, than any other setting 
combination ot these factors. In combinations like high capacities and a movie based 
contract term, the effect of contract term is offset by the relatively loose capacity 
constraint. Similarly, with the 90%/10% contract term, the effect of low capacity 
constraint is offset. 

4.7 Managerial Implications, Conclusion and Directions for 
Future Research 

The proposed GA-based model outperforms the non-capacity based heuristic by 25% 
on an average. This suggests the importance of considering the screen capacities in 
the exhibitor’s decision making. All the main factors, namely, the capacity levels, the 
type of contract tenn, and the number of rapidly decaying movies, are significant. The 
research suggests some important policies for the exhibitors as well as the distributors 
of movies. First, it is implied that the exhibitors should not allocate screens to movies 
based on their demand alone. Even if the movies with high opening strength has 
higher demand in the beginning, its allocation to the high capacity screens might lead 
to inferior schedules because of the contract terms specific to movie types. The 
interaction of factors of capacity and contract term is significant, This suggests that 
the exhibitors with low capacity screens should not blindly exhibit movies of Type I. 
They must take into consideration other factors such as contract temi. 

On the other hand, distributors of Type I movies should go for 90%/10% 
contract term with the exhibitors having low capacity screens while a movie based 
contract term with exhibitors having high capacities. Also the exhibitors with higher 
capacity screens should exhibit Type I movies only at screens with capacities more 
than the demand of other movies if the contract terms are movie based. Exhibitors 
with low capacity screens should prefer movie based contract terms with the 
distributors distributing less number of T)/pe I movies and should prefer 90%/10% 
contract terms with distributors distributing more number of Type I movies. 
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In the current research we liave considered only one theater. However this 
model can be extended to consider an exhibitor who owns multiple theaters. In fact 
this is the common practice, and in any regional area, typically most theater screens 
are controlled by a few chains. Such a setting may allow the theater owner to optimize 
over multiple sites as the decision of scheduling movies on a theaters are highly 
influenced by the other theater. 

The current research can also be extended to address the competitive issues. 
We consider the movies as independent entities and the performance of one movie is 
not affecting the demand of other movie. However in the real world situation the 
performance of one movie may affect the others. Consideration of such effects would 
provide interesting future research ideas. 
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CllAPIERS 

MODEL EXTENSIONS IN 
OTHER APPLICATION AREAS 

In this chapter we discuss the extensions of SilverScreener in addressing some 
problems in similar areas such as television programs scheduling and retail space 
allocations. We also discuss some of the related areas in which meta heuristics such as 
GA can used. 

5.1 1'elcvisioii ProRniui Scheduling 

In their research, Reddy. Aronson and Slam (19^)8) propose a model S.P.O.T. 
(Scheduling Programs Optimally for Television) which explains a network analogy 
with scheduling of programs on television. The model considers an interesting issue 
of 'iend-in” s of television programs. They define lead-in of a program as its strength 
to boost the ratings of a newly introduced or a weaker program. With a minor change 
in model. SilverScreener can also take into consideration this “lead-in" effect and can 
be elTiciently used in television programs scheduling, w ith an advantage of simplicity 
over S.P.O.T. We now explain a model, which is based on SilverScreener and is 
equi\alent to the S.P.O.T. model. We then extend this model for addressing some 
issues, which are not addressed by S.P.O.T. 

S.P.O.T. equivalent model 

We first define the terms to be used in the model; 

N ; Total number of television programs available for scheduling. 

W : Number of uniform time slots to be scheduled. 

LI k ; lead-in effect on next scheduled program due to program k. 

Ij length of program J in number of uniform time slots. 

Tjkw ; binary decision variable 



V’/A» = 1 it' program./ follows program k, and begins at slot w. 

0 otherwise. 

.\7 = binary decision variable 

Xj,w =1 if program 7 is scheduled on slot vv. 

Pj ,, contribution to objective function if program j is scheduled at slot w. 


The S.P.O.'f. equivalent model will be as follows: 


N ir 
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C'onstraint (5.1) checks that only one show is scheduled at any point in time. 
Constraint (5.2) takes into account the lead-in effect of show k, only if it is scheduled 
just before show /. 

To illustrate the model developed, we explain a simple example. Suppose that 
there are 10 programs to be scheduled and we have a total of 6 slots. Hach program 
has the play-Icnglh of 2 slots. Program number 7 has a lead-in eifect and is scheduled 
in slot number 2. "fhus now if we have to schedule program 3 at slot 1 then 

'I'hc constraint (5.1 ) will give for slot 2, 

X^j + .Y ^ j 4- 4" .Y 7/ 4* .Y 72 1 





i e- ()+()+l<= 1. which is not satisfied. Thus Constraint (5.1) will limit 

the number of programs to he scheduled at any slot to 1 . 

Now if program 3 is scheduled at slot 5, we have .\-.{. 5 =l. Then by constraint (5.2), 
(5.4) and (5.6), we have Vi7.5 < = {x 3 s + X 73 )/2 i.e. <= ( l+0)/2 = 0. Thus as program 3 
did not immediately follow the program 7, the lead-in effect of program 7 will not be 
considered. This way both the constraints work well and the above model is 
equivalent to the S.P.O.T. model. 

We can add one more practical consideration into the model. The S.P.O.T. 
model does not incorjxirate daily shows. We present a model that considers the 
telecasting of programs lliat are shown daily throughout the week. We define 

/%,.• Contribution of show / towards objective function, when scheduled at time slot 
w of day /'. 

.V;„, .• Binary decision variable takes value 1, if show ; is scheduled at time slot w of 
day /., 0 otherwise. 

Bjw : Contribution of show / towards objective function, when scheduled at time slot 
w for the whole week. 

Yjkwi : Binary variable, takes value 1, if show j is scheduled at time slot w on day i, 
and show k is scheduled just before that, 0 otherwise. 

Zjw : Binary variable, takes value 1, if show j is scheduled at time slot w for the 
whole week, 0 othenvise. 


The model will be; 

N H' 7 


N N ly 7 


N ir 


/=i w=i i=-i y=i *=! w=i /=! y=! H'=i 


f» ■» 


Subject to 

y Jhvi - jk»i "*'■ ^ jw '^J,k,W,i (5.7) 
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Constraint (5.7) ensures that lead-in effect is considered only when showy follows 
show k. Constraint (5.8) ensures that at one time slot on a particular day only one 
show is scheduled. Constraint (5.9) insures that over a week at a particular time slot 
only maxiinuni of 7 shows are scheduled. 

We illustrate the above model with an example. Suppose there are 10 slots per 
day to be scheduled and there are 30 programs with a play length of 3 slots each. If 
program 7 is scheduled for a week on slot 6 then Z 7 fi=\. Now if we schedule program 
3 on slot 7 of the day 4 then .v.!?./=l . 

Constraint (5.8) gives 

.\\r 4 + .\yrtj...+ .v.iw *- 7(5 < = 0. i.e. l+....0.+l< = 0, which is not 

satisfied. 

Also by constraint (5.7) the lead-in effect of program 7 will come into effect 
for program 3 only when either: a) program 7 is scheduled exactly before program 3 
on one day, or b) program 7 is scheduled as a daily program and program 3 follows 
program 7 on one day, or c) program 3 is scheduled as a daily show and program 7 
precedes program 3 on one day, or d) both programs are scheduled as daily shows and 
program 7 precedes program 3 on all 7 days. 

This model still does not consider a show’s re-telecasting possibility. If we 
consider re-telecasting, the model will be nonlinear as each re-telecasting will affect 


the contribution of show towards the objective function. This extension will thus 
require a heuristic method such as the one explained in previous chapter. The 
proposed Genetic Algorithm heuristic can directly be used for solving this complex 
problem since with the possibility of re-telecasting, it can also consider the 
deterioration of contribution to the objective function value with each re-telecast. 

5.2 Retail Space Management Considering Demand Elasticities 

Corstjens and Doyle (1983) explain a model for retail space management. Their 
model considers elasticity in demand (sales) of a product with varying amount of 
space allocation to it. They also consider a cross elasticity in demand of products. 
That is, they consider the eriect of certain space allocation of one product to the sales 
of another product. Such conditions are interesting and GA can be applied to such 
problems, which will help obtain near optimal solution. For the similar problem of 
retail space management, Boris, Farris, and Freeland (1994) explained the use of 
another meta heuristic. Simulated Annealing, after modelling the problem as a 
constrained optimization problem. 

5.3 Internet Space Scheduling 

The number of World Wide Web users is increasing at a rapid rate. This has made 
Internet an important media for advertisement. The scheduling of adds on the web site 
poses some interesting problems in front of web site administrators. The pricing is 
based on the site exposure or the banner ad exposure. An advertisement on the web 
has two parameters, namely the size of advertisement and its display frequency. To 
schedule ads on the web in a way as to optimize on both the parameters is a difficult 
task. Kumar, Jacob, and Srikandarajah (2000) use GA in this area. They propose an 
algorithm based on GA, which they call hybrid-GA that gives significant increase in 
objective function values as compared to other existing algorithms for this problem. 

This problem of advertisement scheduling is analogous to the retail space 
allocation problem discussed earlier. The site administrator has certain space, which is 
to be allocated with an objective of maximizing the revenue earned. At each decision 
epoch, the administrator takes a decision of which advertisements to schedule and 


what space to be allocated to each advertisement. Depending upon the action taken by 
the administrator, the system makes transition to a new state. This state can be defined 
as the function of the advertisements scheduled and the respective space allocations. 
A piobability distribution of revenue that will be generated is associated with each 
state. Thus the Internet ad scheduling problem is in one way a combination of the two 
retail space allocation problems we have considered, which has a complexity in terms 
of varying space allocation to each advertisement as well as stochasticity of revenue 
corresponding to each allocation. 


CHAPTER 6 


CONCLUSION AND LIMITATIONS 

In this thesis we applied marketing and management science (operations research) 
methods along with a tool of Genetic Algorithm to dynamic retail management of 
movies. We showed how quantitative model building and meta-heuristics like Genetic 
Algorithms can assist marketing decision making in complex environments involving 
perishable products. We addressed the general retail space allocation problem, in 
complex situations such as stochasticily in the demand of products and the display of 
multiple products. This space allocation problem was studied in the context of motion 
picture exhibitors. 

We described the exhibitor problem in two ways. In the first problem the 
exhibitor has a single screen and the demand distribution for movies is considered to 
be stochastic. We propose an MDP based model to address this problem. In the 
second problem we consider an exhibitor with multiple screens and assume 
deterministic demand. We propose a GA-based heuristic to handle these problems. 
We explain how the models might assist the movie exhibitors in making profitable 
scheduling decisions. At a broad level, the model would benefit managers in two 
ways. One is by helping managers improve their scheduling techniques. That is, the 
model could help the manager adopt a new decision making style which is based on 
an analytical approach or a tested heuristic. The second is by helping the managers 
perform the same scheduling task in a more efficient way, using computer-based 
routines. 

Although we have considered movies in this research, our results can readily 
be applied to any other perishable entertainment product. Thus in a broad sense, the 
current research addresses the problem of general retail space management. We also 
present extensions of the SilverScreener for its application in some other similar 
areas. We propose GA-based heuristics for addressing more complex situations, 
which cannot be efficiently tackled by simple models. 



The current research has some limitations that suggest attractive ideas for 
future research. We address a single-screen problem with stochasticity in demand, and 
also a multiple-screen problem considering deterministic demand pattern for the 
movies. The problem of multiple-screens considering stochasticity in demand which 
will be very close to a practical situation, is yet unanswered. We believe that this 
problem can be answered efficiently a combinations of heuristic methods including 
GA. 

There are computational limitations to our problem. The MDP code developed 
has limitations for the size of problems handled. We used certain state reduction 
techniques to solve the problem of sixteen movies in a scenario. The time take by this 
code for a 16-movic problem is about 45 minutes, which is very large. For the 
application of this model in the real world, it should be able to handle larger 
problems. A user friendly Decision support system based on these models can be 
developed. 

To summarize, we present application of MDP approach in the area of retail 
space management, and the use of GA approach in the field of marketing. The thesis 
attempts to bring together the two important fields of marketing and management 
science/operations management. 



APPENDIX 1 


MDP CODE IN C LANGAUGE 


// This program uses 3 ratiks and for 16 movies 
single link list 
#include<stdio . h> 

#include "prmtrsoml6 . h" 

// defining structure for a state 
struct state 

{ 

int stag; 

int indl; 

int lev; 

int rumovrank; 

float maxm; 

int max_ind; 

struct state *fptr; 


this program uses 


void readataO; // function which reads data 
void mat_gen ( ) ; 

//generated matrix of 3x8 for all combinations of 1,2,3 for 8 movies. 

void consetO;// generates consideration set foreach stage. 

int check_state_f easibility (int mov,int lev) ; 

float get_prb(int l,int levl,int stage, int action, int runmv,int 

act, int movl); // gives probability of transition from one state to 

other 

ma i n ( ) 


struct state *list , ’^'listl, ^current; 
int 

stage , i , mov, run, levl , j , act , c, 1 , action, t ,posl, m, k, runmv, flag, count , 
reduce ; 

float MAXIM, tran_prb, X, exp=0 ; 

FILE ‘^fout; 

f out = f open ( '* resconB . dat " , '^w" ) ; //output file. 
list= (struct state *) malloc (sizeof (struct state) ); //this will 
be the header for link list 

readata 0; //reading data from file 

mat_gen();// generating combinations of ranks from 1 to 4 and 
storing in array level [NUM_LEVELS] [NUM_M0VIES] 

conset();// generating a consideration set for each stage and 
storing in conset [NUM_WEEKS] [4] 

for (i = l ; i<^NUM_M0VIES/2 ;i + + ) con_set [0] [i] -con_set [1] [i] ; 
current =1 i St ; //current record is initially the header; 
for (stage«NUM_WEEKS; stage>=«l ; stage--) 

//loop for all currently running movies which are in the 
consideration set of stage-l; 
for (i=X ; i<xNUM_M0VIES/2 ; i++) 

^ mov=con_set [stage- 1] [i] ;// considering a current movie, which 
will be one of the movies in action set for previous wee . 
if (movl «0) 


{ 

;run..)//for each possible 

{ 

for (k = l;k< = 3;k+ + ) //for each possible rank 

f lag=l ;posl=0; 
for (m-l;m<=NUM_M0VlES/2;m++) 
if (con_set[staie] [m] ==mov) ( flag=0;posl=tn- } 
for{l = l;l<=NUM_LEVELS;l + + ) i 3 -F n, / 


if ( (flag= = 0) && (level [1] [posl] !=k) ) continue; 
if(k!=4) //checks the feasibility of the state 

current ->lev=l; current- >indl=mov* 10 +run;current->maxm =0 ; 

current - >max_ind=3 0 ; current - >stag=stage ; current - >rutnovrank=k • 
for (act=l ;act<=NUM MOVIES/2 ;act++) 

{ 

MAXIM=0; 


action=con_set [stage] [act] , -//movie 
if (actionl “0) 


{ 


to be considered as action 


c«level[l] [act];// current rank for movie selected as 

action. 

if (act ion«amov) t=run; // if action is currently playing movie 
else taO; 

MAXIM+arev[c] [t] ; 
if (stageUNUM^WEEKS) 

( // freeing records for stage next to next since they are 
no more required. 

if (stage<NUM_WEEKS - 1) 
while (list->stag>stage-i'l) 

( 

listl-list-'>fptr; 

f ret;‘ { 1 i liit ) ; 

J i, li t ~ I i t;5 1 1 ; 

} 

listl=:list; 

// moving state by state in records of next stage to calculate 
revenue* transition probability, which will sum up to revenue 
corresponding to this action. 

while (listl->stag! =stage) 

{ 

if ( list 1 - >indl-=action*10+t + l) 

{ 

1 ev 1 . 1 i ut 1 ■ : > h .‘V; 
runmv=listl“ >rumovrank; 

tran_prb=get_prb ( 1 , levl , stage , act ion, runmv, act , mov) 
//getting probability of transition from current level to next level 
MAXIM+=tran_prb*listl->maxm;// added to total revenue 
corresponding to this action. 

} 

listl=listl->fptr; 

} 

} 


if (MAXIM>current - >maxm) 

{ 

c 1 1 r r e n t “• > n^a x m === M AX I M ; 
current - >max_ind:=act ion / 


} //printing out current stage, movie playing, run length, rank of all 
movies, and optimal action and optimal revenue, 
if ( (stage=.^l) && (mov= = l) ) if (current->maxm > 0 ) 

{ 

printf ("\nstage %d movie %d run %d rank %d ranks" , stage, mov run k)- 
f or ( j =1 ; j <=NUM_M0VIES/2 ; j + + ) if {con_set [stage] [ j ] ! =0) printf ( “ Id" ' level 
[1] [j] ) ; 

printf ("maxmoney If optact Id" , current- >maxm, current->max ind) ; 
x=current “ >maxm; “ 


for ( j =1 ; j <=NUM_M0VIES/2 ; j + + ) if (con_set [stage] [j] [ =0) x=x*initial pr [co 
n_set [stage] [ j ] ] [level [1] [j] ] ; 
exp=exp4 X; 
f print f ( font , " ) ; 

} 

if (current - >maxm>0) 

{ 

current->fptr= (struct state *) malloc (sizeof (struct state)); 
current ^current - >f pt r ; 

} 

} 


count =0 ; reduces 1 ; 

for ( j “ 1 ; j <sNUM_M0VIES/2 ; j + + ) { if (con_set [stage] [j] ==0) count=count + l ; } 
for { j =1 ; j <scount ; j + + ) reduce* =4 ; 

Isl+reduce- 1 ; 

} 

} 

} 

} 

} 

} 

fprintf (fout, "EXPECTED REVENUE lf",exp); 
return (0) ; 

}//end of main; 


//function for reading data from file 


void 

( 


readata ( ) 


FILE *infile 
int j , k, 1 ; 
if ( (infile * 


fopenC’dataslS.dat”, "r”) ) == NULL ) 


pr int f (" could not open 


datal.dat ' ") ; 


exit ( 1) ; 


}; 

/* the immediate rewards for any movie */ 
for(j =. 1 ; j <= NUM RANKS ; i+ + ) 

( 

for(k = 0 ; k < NUM_WEEKS ;k++) 

fscanf (infile, "%f'',&rev[j] [k] ) ; 

^ piintfC'j %d k %d %f ",j,k,rev[j] [k] ) 

printf ("\n") ; 

} 

/* assigning initial probabilities */ 
for(j = 1 ; j <= NUM_MOVIES ; j++) 
for{k = 1 ; k <= NUM RANKS ;k++) 

{ 

tucdiit (i utile, "%t" ,&initial_pr [j] [k] ) ; 


/* assigning transition probability */ 


for(j «-• 1; j <= NUM_MOVIES ;j + + ) 

for{k = 1 ; k <:= NUM_RANKS;k+ + ) 

ford = 1; 1 <= NUM_RANKS ; 1 + + ) 

fscanf (infile, "%f " , &ind_pr [j] [k] [1] ) ; 


/* release dates */ 

for (j ^ 1 ; j NUM MOVIES ; j+ + ) 

( 

fscanf (infile, "%d" , &rel_dt [j] ) ; 

} /’^printf (" rel dates %d, %d" , j , rel_dt [ j ] ) ; } */ 

/* obligation periods */ 

for(k ^ 1 / k NUM MOVIES /k++) 

{ 

fscanf (inf ile, "%d" , &OPD [k] ) ; 
printf ("opd, %d\n’' ,0PD[k3 ) ; 

); 


/* TYPE OF MOVIE*/ 
for(k=l;k<= NUM_MOVIES; k++) 

( 

fscanf (infile, nd",&TyPE[k] ) ; 

} 

f close (infile) ; 

} 

///function to return probability of transition between states, 
float get_prb{int 1 , int levl,int stage, int action, int runmv, int 
act,int movl) 

( 

float prob*l; 

int mov, prev_rank, new^rank, pos, f lag2; 

int i, j, flag; 


for ( j.--=l ; i<=NUM__M 0 VlES /2 ;ii- + ) 

{ 

, B,.t- [stage! [i];// using "stage" variable to get 

con s i c;lc l a 1. 1 on b t j t , a nd 90 i ng mov i e by mov i e . 

1 f ( mov.» 0 ) cont. i riue ; 

^ prev^rank = level [1] [i] ;// getting previous rank of movie using 

if ( (mov! -movl) &Sc{mov==action) && (prev_rank==3) ) prob=0 • 

pOS^^O; ' 

flag2=0; 

if (mov3«movl) flag2=l; 

for ( j=:l; j<=NUM_M0VIES/2; j+ + ) 

{ 

if (con_set [stage + 1] [j]==mov) 

{//finding position of movie in next stage's consideration set 

pos = j ; 

break; 

} 

} 

if (pos ! 3^0) 

{ 

new__rank=level [levl] [pos];/ getting new rank, 
if ( (flag2 = l) &&(mov!=action)&&(new_rankl=3) )prob= 0 ; 

//if movie replaced and rank still not 3, prob=0. 

else if ( (flag2 = l)&&(mov!=action)&&(new_rank==3))prob=prob; 
//if movieis replaced and rank =3, then prob=l. 

else prob*3ind_pr [mov] [prev_rank] [new_rank] ; 

//else use of individual transition probability. 

} 

} 

// this was for the current running movie and action, now for other 
movies . 

f or ( i = 1 / i < =NUM_MOVI ES/ 2 ; i + + ) 

( 

flag-1 ; 

mov-con_set [stage+1] [i] ;//one movie at next stage is 

taken 

if (mov===0) continue; 
prev__rank=level [1] [i] ; 
pos=0 ; 

for ( j - 1 ; j <=;NUM___MOVIES/2 ; j ++ ) 

{ 


i f (con_set [stage] [j]==mov) 

{ 

f lag=0 ; 
break; 


}//if movie newly entered consideration set, its 
initial probability is considered, 
if (flag=.‘-l) 

new rank^level [levl] [i] ; 


prob*=initial_pr [mov] [new_rank] ; 


} 

} 

flag 1 ; 


for ( i ^ - 1 ; i< NUM MOV I 

i f { f 1 ag 1 ) proh 

retin n ([nol)) ; 


f (con set [stage.u [i] ==action) flag=0; 
ind p, (action] [level [1] [act]] [runmv] ; 


} 


} 

//function generating matrix of rank combinitions . 
void mat__gen (void) 

{ 


int q==l; 
int i , j ; 

int mO, ml, m2, m3 ,m4, m5, m6,m7; 


for (mO=l ;mO<=NUM_RANKS;mO++) 

for (ml = l ;ml<=NUM_RANKS;ml++) 

for (m2=.l;m2<=NUM_RANKS;m2 + + ) 

for (m3=l;m3<=NUM_RANKS;m3++) 

f or ( m4 « 1 ; m4 < »NUM_RANKS ; m4 + + ) 
for (in5 = l ;m5<=NUM_RANKS;m5++) 

for(m6=l;m6< =NUM_RANKS ; m6 + + ) 
for (m7=l;m7<=NUM RANKS, •m7++) 

{ 

level [q] [1] =m0; 
level [q] [2} =ml; 
level [q] [3] =m2; 
level [q] [4] =m3 ; 
level [q] [5] =m4 ; 
level [q] [6] =m5 ; 
level [q] [7] a^m6/ 
level [q] [8] =m7 ; 
q++ ; 


// function generating the consideration set at each stage using the 
data of released movies, movie types since type 3and 4 movies do not 
stay in consideration set after 4 weeks, 
void consetO 

( 

int i , j , k, goodmv, stage, flag, 1, flagl; 
int topset [NUM_MOVIES+10] ; 
for (i»0 ; i<«NUM_MOVIES; i + + ) topset [i] »0; 
j »1 ; 

f or ( i«l ; i<=s NUM MOVIES ;i + + ) 

{ 

if (TYPE [i] -»1) 

topset [j ] “i ; //storing alltypel and 2 movies. 



3 > + ; 


} 

f or ( i = j ; i < -NUM_MOVIES ; i + + ) topset [i]=0; 

for (3tage-NUM._WEEKS,-stage>=l; stage--) //for each stage. 

for { i .^NUM MOVIES ; i >= i ; i - . ) 

{ f:iagl=:l; 

f o r ( j 1 ; j < .-.=NUM_MO VI E S ; j + + ) 

if (topset [j] = = i) 

{ 

f lagl=0 ; 
break; 

} 

} 

i f ( ( f lagl=;=0) ScSc (rel_dt [i] <=stage) ) // if movie is good and 
released then added to consideration set 
{ 

con_Bet [stage] [k] «ii ; 
k t f ; 

} 

}// else writing movies in con set until stage/2 or 8 movies 
are written. 

i f ( ( (8tag€^>^«NUM WEEKS/2) && (k<=:NUM_M0VIES/2) ) 

1 I ((stage<;NUM WEEKS/2) && (k< = stage) ) ) 

i f (Btage<NUM WEEKS/2) j~stage;elae j =NUM__M0VIES/2 ; 
while (k<^ j ) 


f or ( i ES ; i ; i - - ) 

{ flag-l; 

f or (1 ; l<;rk; 1 + + ) 

{ 

if ( { ir^=^con set [stage] [1] ) | | (rel_dt [i] >stage) ) f lag=0; 
} 


if (flag --1) 

( 

Gon_set [stage] [k]«i; 
k4* ^ ; 


} 

}}} 

for ( i»NUM_WEEKS ; i >»1 ; i - - ) 

{ printf ( " \nstage %d",i); 

for ( j>=l; j<=NUM_M0VIES/2; j+ + ) 

{pr intf ( "\t%d" , con_set [i] [ j ] ) ; 

} 


getcharO ; 


} 


} 



in^.ADI^^R li'lLK FOR IMDP MODEL 


#ifndef PRMTUS 
#define PRMTRS 
((define NUM_ WEEKS 8 
((define NUM MOVIES 16 
#define NUm”levelS 6561 /* 
of movies*/ 

# define NUM RANKS 3 


from combination of no. of ranks and no 


int con_set [10] [NUM_MOVIES+10] ; 

float rev[NUM_RANKS + 10] [NUM_WEEKS+10] ; /* revenues in a state and 

week */ 

float initial_pr[NUM_MOVIES+l0] [NUM_RANKS+10] ;/* initial orob */ 

float ind_pr [NUM_MOVIES + 10] [NUM_RANKS+10] [NUM RANKS + 10] •/* 

individual probabilities */ ~ ‘ 

int rel__dt [NUM_MOVIES+10] ; /* release dates */ 
int level [NUM_LEVELS + 100] [NUM_MOVIES + 10];/* matrix to be 
generated*/ ” 

int OPD[NUM_MOVIES + 10] ; /* obligation periods */ 

int TYPE [NUM__MOV1ES + 10] ; 

#endi f 



appendix - 2 

CA-(X)DE IN C++ LANGAUGE 

#inc] ude* Hldi o . li > 

#incliide<Bt d 1 i i) . h > 

#incl udectuat: h Ai> 

#define NUM^MOVIES 38 
# define NUM^^WEEKS B 
#define NUmJsCR 6 
#define MAXSLOT 48 
#include<graphics .h> 

# i nc 1 ude < con i o , ii > 

# i n c 1 u d e < i oina n i p . h > 

if i nc 1 lid (.3 < a t i' i . h > 

# i nc .1 ud(* < f B t i vrim . h > 

//#detine NUM^ U 

//defining cluomosome an a structure with b as an array of schedule 

Struct: tdurom 

( 

int no; 

int b (NIIM WEEKSUNUM^SCR] ; 
f loat, revenue; 
i nl copies ; 

St: ruct chrom *next ; 

); 

void sel ect: ion ( ) ; 
void c r os sever (); 
void mutatiouO ; 

float eva 1 uat ion ( int b [NUM^WEEKS] [NUM^SCR] ) / 
int 

rel^dt [NUM^MOVIESl .demand (14] [NUM_WEEKS+10] .movie^type [NUM_^MOVIES+10] 

f 

float contract_term[14] [NUM_WEEKS+10] ; 

int rrr. acreen^cap [NUM_5jCR]7 

float a [NUM^MOVIES] [NUM^WEEKS] [NUM^SCR] ; 

struct chrom ^head, *h; 

int popsize . maxgen; 

float pc . pm. tot_popO , tot_popl . tot_pcO . tot_pcl . tot_pmO , tot_pml ; 
int gen-0; 
mainO 
{ 

FILE *fp; 

int ab, be, cd, de, hpop»0,best_pop; 

float higheat^O . 0 , hpe^O , hpm«0, tot_pop, tot_pc, tot_pm, best_pc, best_pm 

static int beet^revenue [ 3 ] [ 3 ] ; 

int gdr i ver»DETicT, gmode*0 , X, y, z , k, xx, count ; 

int n, 1 , f lag, f lagl , f lag2 , r ,m; 

struct chrom *p,*q,*best; 

f pas f open ( »'scenl6 .dat *’ , "r” ) ; 

initgraph (agdriver , &gmode, " ) ; 


//giving screen capacities 



screc^n cap [0] »3472 ; b 
] --12 08 / Bcieun ccdp ( 
coni ract i (m iu 1 1 ) ( 1 ) 

. 1 ; 

con t r a c t . t; ( * i/m [ 2 1 ( 2 1 
. 4 ; conLr\u.:t. tern\ [4 ) 
for (x . 1 ; ;Xi f ) 


1 oen cdp [1] =2736;screen_cap[2] =1728 /screen cap [3 
> I 1 1 12 ;screen_cap [5] =904 ■ ^ieen_capi3 

). l;contract_ter.n[l] [2] =o:3;contract_term[2] [l]=o 

).2;contract_term[3] [1] -0 . 2 ; contract term[3] [2]=0 
1| -0.2;contract_term[4] [2]=o.4; ~ l j u 


for (y- 3 ;y<-.:NUM WEKKS/yt + ) 


if (xa = l) contract_term[x] [y]= 0 . 5 ; 
else if (x== = 2) contract_term[x] [y]= 0 . 3 ; 

else contract_term [x] [y] =0 5; 

} 

for (x=0;X<NUM MOVIES;x++) 

( 

fscanf (fp, "%d",&movie_type[xl ) ; printf (■' 

% d , mov i e t ype ( x ] ) ; 


{vrint f: { ) /tjcl rh {) / 


//writing demand for type of movies in each week 

demandtl] [0] ; demand [ 1] [ X] ^-2189 /demand [1] [2] =1306 ;demand[l] [3]=78 

0; 

demand [ 1 1(4] » 4 G ; dettui nd 1 1 1 [ S ] 277; demand [ 1 ] [ 6 ] = 1 6 5 ; demand 1 1 ] [ 7 ] = 9 9 ; 
demandtal (0) 41 02 /demand [2] [lba280;demand[2] [2] =2622 /demand [2] [3]l20 
97; 

demand(2] ( 4 1 ^.1677 ; demand [2] IB] «1341 ;demand[2] (6] =1072 ;d€mand[2] [7]»85 
7 ; . 

demandO] [01 -1804 /demand [3) [ 1 ] =1157 /demand [3] [2] =742 /demand [3] [3] =476 
/ demand[3l [4] «30B/deniand[3] [6] =196 /demand [3] [6] =126 /demand [3] [7]=80; 
demandl4l to] :«1S77 /demand [4] [ 1] =1358 /demand [4] [2] =1169 /demand [4] (3]=10 
07; 

dernandMl (4) -867 ;demandl4) 15) .»747;demand(4l [6] =643 ;demand[4] [7]=554; 
for ( 0 ; xcNUM MOV 1 KS ; X f f ) { 

fscanf ( fp," %d",&rel_dt [x] ) ; 
printf ("%d " , rel_dt [x] ) ; 

}getch ( ) ; 

for (ab-0 ; ab<2 ; abf + ) 
f o r ( be 0 / be < 2 ; be f f ) 
for (cd--0 / cd<2 ; cdf t ) 

foi (de*d);de<4;deM ) //for loops for different settings of 
parameters 


gen-O; 

haad« (struct chrom *) malloc (sizeof (struct chrom) ) ; 
ti- liead; cl i ( ) / 
switch (ab) 

case 0 : popsize=3 00 /break; 
case 1 ;popsi2e=400;break; 

} 

switch (be) 

case 0 : pc» 0 . 9$ /break; 
case 1 : pc* 0 . 99 /break; 



Bwi, t. cll (ccl) 

{ 

case 0:pm=0.03;break; 

^ case 1 :pm=0.05;break; 

Bwi tdi (de) 

{ 

case 0 :rrr=l; break; 
case 1 :rrr=25, -break; 
case 2;rrr=100;break; 
case 3 :rrr=400;break; 

^ case 4;rrr=6000;break; 

maxgen=2000 ; 
pc*=100; 

gotoxy{10,10) ;printf ("give random seed ") .-//getting random 

seed for random number generator, 
scanf ( ’'%d” , Scrrr) ;clrscr () ; 

Brand (rrr); 
pm*«100; 

//creating initial population 
f or (x«0 ; x<popBi ze ;xi f ) 

{ 

h*“>next» (struct chrom ’^) malloc (sizeof (struct chrom)); 
if (h- >next»i^sNULL) (printf ( "could not assign" ); exit ( 1) ; } 
h“ >next ->noax/h'* >next">copies=0;h->next->next=NULL; 
for (y=«0;y<NUM_WEEKS;y+ + ) 
for (Z!a0 ; z<NUM_SCR; Z + + ) 

( flag®0; 

while ( f lag»*0) 

{ 

flagl-O; 

r=rand()% NUM_MOVIES; 

for (1 = 0 ; l<z if (r==h->next->b [y] [1] ) flagl=l ; 

//considers release date for initial population generation, 
if ( (rel_dt [r] <=y) &:&(flagl==0) ) 

{ 

h- >next- >b [y] [z]«r; 
f lag=l ; 

} 

} 

///using repairing function to make chromosomes feasible, 
for (xx»0 ;xx<NUM_SCR-2 ;xx+ + ) 
for (y*2;y<NUM_WEEKS;y++) 
f or (zwO ; z<NUM_SCR; Z++) 

{ flagl«l; 

rah->next“>b [y] [z] ; 
f 1 ag « 1 ; 

f or (k»0 ; k<N'UM_SCR; k+ + ) 

if (r»wh" >next->b [y-1] [k] ) (f lag= 0 ; f lagl=0 ; } 
if(fla 9 ««l) 

{ 

f or (m»0 ;m<y“ 1 ;m++) 


f or ( k 0 ; k < NUM_SCR ; k+ + ) 


^ i t { i ^ >next- >b [m] [k] ){ flag=o , -break, • } 

// chc;t:k i ng tor di acoiit inui ty of movie in b[yl [z] • 

1 f. ( (f lag- - 0) && (f lagl = = i) ) //movie discontinuous 
{ count- 0; 

for (m-^0;(iKy-l;m-f + ) 
f or ( k == 0 ; k<NUM__SCR ; k-h-f ) 
if (rr^=h->next->b[m] [k] ) (count-H-f ; } 
k-0; 

f lag2=0 ; 

if (count>NUM_WEEKS/2) 
while (f lag2==0) 

{ 

f lagl-0; 

n-^h->next->b [y-1] [k]; 

for (1 :0;1<NUM_SCR;1++) 

if ([i h .•m;xt >b[yl[l])flagl.l; 

it (1. lai),1....0) 

{ 

h - >next->b [y-l] [k] =h->next->b[y] [z] ; 
f l aq2-l ; 

1 

Km ;// BchtHiuling movie in the weeks where it was not scheduled 
r emov i rit| d i ecun t: i nu i t y . 


else 

while ( f lag25SKO) 

{ //rc^moving the movie, because it was for less number of times. 

f lagl»0 / 

n-h»>next«>foty“l] [k] ; 
tor ( 1 ^ 0 ; 1 <]MUM_SCR ; U O 
i f >next « >b [y] [11 ) f lagl=l ; 

i t { f 1 ag I ^ 0 ) 

{ 

h- >next- >b [y] [ 2 ]«n; 
f lag2 » 1 ; 

} 

k-H-f ; 


h“ >next '• >revenue»evaluat ion (h->next->b) ; //evaluating chromosome. 

h-h“ >next ; 

) 

hahead; 

h»head;best» ( {struct chrom *) malloc (sizeof (struct chrom) ) ) ; 
best - >revenues«h» >next - >revenue ; 
for (yaO;y<NUM_WEEKSfy++) 

for (2«0 ;z<NUm"”sCR;Z4-«*-) ^ ^ , 

besC->b[yT[z]«h->next->b[y] [z] ;// finding the best schedule 

for (x»0 ;x<popBize;x-t--f) 

if (h->next->revenue>best->revenue) 


r 

I 




I 


best - >revenue=h- >next->revenue • 
for {y::.0;y<NUM_WEEKS;y++) 
lor (h:=^0;2;<NIM_SCR;Z + + ) 

^ h(M)t, >blyl [z] =h->next->b[y] [z] . 

h h >iiext: ; 

} qol oxy { 10, 10 ) ; 

print!: ("\nb«iit revenue .. " ,beBt->revenue) • 

h - h e a d ; c 1 1 ' a c r ( ) ; ' 

gotoxy (20,20) ; 

printf('PltBaB 0 wait, till I di splay ths ire suit ” ) * 
for (k^0;k<maxgen;kf+) ' 

{ //moving generation to generation 

srand(rrr); 

selection ( ) ; 

Brand (rrr) ; 
croeeover ( ) ; 

Brand (rrr) ; 
mutationO ; 
h-l'iead ; 

p 4iead >next ; q* head - >uext ; 
tot (x-0 ; xcpopsi ze ;X+- + ) 

{ 

if (h >next->revenue>p->revenue) {p=h- 
>nex t / } / / f i nd i ng l)ea l 

it (h >next“>revenue<q->revenue) {q=h- 
>next ; } / / i i nd i ng worst 

hish™ >next ; 

} 

if (beat - >revenue>q" >revenue) 

{ 

q->revenue«be8t->revenue; //replacing worst by previous best, 
for (y-0;y<NUM_WEEKS;y++) 
for(z»i0;z< NUM_.SCR ; z + + ) 
q- »b (yl [z] *:best- >b [y] [z] ; 
if {best - >revenue>p->revenue) 

{ 

p- >revenue«‘beat - > revenue ; 

£ o r ( y - 0 ; y < N UM,_WE E KS ; y + + ) 
f or ( Z == 0 ; z < NUM~SCR ; 2 + + ) 
p - >b [y] [z] «best->b [y] [z] 

} 

else ( //replacing previous best for current best 
best - >revenue»p- >revenue ; 
for (yaO ;y<NUM_WEEKS;y+ + ) 
f or ( z « 0 ; 2 < NIJM^SCR ; z + + ) 
best ->bty) [2]«P'*>b[y] [z] ; 


if ( ( (k+l) %100-«0) j 1 (k==maxgen-X) ) 
clrscr 0 , -gotoxy (15/ 15) ; 

printf("%d %f %f",popsize, pc,pm) ;gotoxy(l5, 15) 

printfC'best revenue » %f , generation %d",best- 
>revenue, k+l) ; 

} 


i t ( (k itiaxtjuii 1 ) ) 

( 

[11 i lit, I ( " \m" ) ; 

iui (y ();y' NUM Wlih;KS;yn ) 

{ t '» (■■■ U;i'.<NUM SCR;2: + f) 

id", beBt->b[yl[2l)- 

printt:("\n") ; 


gerif' 4- ; 

i f (best ’• >revenue>highest) 

( 

highest =best“ 

>revenue;hpop=.popsi7.e;hpc:.pc;hpm=pm;//changing highest for current 

setting. 

} 

} 

bent^revonue [ab] [be] +abest->revenue; 

P'.head;q head; 
whi le (q! NOI.b) 

( 

p-q- >next ; 
free (q) ; 
q^p ; 

} 

free (beat) ; 

} 

tot pop 0«0 ; tot, popl«<0; tot__pc 0 = 0 ;tot_pcl= 0 ;tot_pm 0 = 0 ;tot_pml= 0 ,- 

f or ( t)C» 0 ; bc< 2 ; bc+ + ) 
tor (cd^iO ; cd<2 ; edf + ) 

{ 

tot_pop0 + «be8t_revenue [0] [be] ; 
tot popl4-ssbest_revenue [1] [be]; 

1 

for (ab»«0 ;ab<2 ;ab+f) 
f or ( cd a 0 ; cd< 2 ; cd+ + ) 


tot_pcO + abest_revenue [ab] [0] ; 
tot_pcl4-«best_revenue [ab] [1] ; 

} 

for (ab^O ;ab<2 ;abf +) 
for (bc-O ;bc<2 ;bC4*4*) 

tot_pmO + «best_revenue [ab] [be] [0] ; 
tot jml'Habest^revenue [ab] [be] [1] ; 

if (tOt^;pop0>tOt_popl)best_popa300;else best_pop=400; 
if (totlpc0>tot3c:l)best_pca0.95;else best_,pca0 . 99; 
if (tot_pmO>tot_pml) best^pmaO . 03 ;else best_pma0.05; 

printf ( “\n’‘ ) ; 

printf('’ best popsize«%d/ bestpc=%f, 
be8tpm«%f , best pop, beBt_pC;best_pm) ; -ri u- 

print f { “ Xnhighest reva%f , hpopa%d, hpc=%f , hpma%f , highest , 

hpop, hpc, hpm) ; 

closegraphO ; 
geatchO; 


} 


i Vi Mill ( 0 ) ; 


void 1 (mM: i t »n ( ) //:;<' i vci i on process 

( 

int M<‘ 1 iMAXiJI.uM' I I HOC) j , 1 I at| , eelection iMAXSLOT+isoo] ‘ 

int X, y . -w I , II, k ; 
struct: t:hi nm ‘^p, *q; 
p»head * >iu‘xi. ; 
for (n«d);n<2;n4* + ) 

{ 

sapope i ze; p-head- >next / 
x-^0; 

while (pi *NULL) 

{ 

p- >no-^x ; 

X f ♦ ; 

■ pr-p ;>nexl ; 

} 

i «r ( X 0 ; x< pcqui i ze ; x m ) eel [x] wx; 
f or { K - 0 ; x-< popu i ze 2 ; x f + ) 

( 

I • raud ( ) % z; 

Uidi'Ct ion [xl -oKil [rl ; //chromosomes selected randomly 
i vn (y r ;y*;z ;y + 1 ) ael [y] asel [y+l] ; 

Z ■■ - 1 ; 

} 

select u>ulxl [01 ; 
ae I ect i ou [ x 4 1 1 « ae 1 [ I ] ; 
f u i: { X 0 ; K < P<1| J u i z e 1 ; X r r ) 

{ y select ion (x] ; 

p htsid;qd3iead; 

wlii ]e(p 5*next->nol «y) {pwp->next;} 
ut‘ I ection [x+l] ; 

whi le {q* >next->nol ^z) (q«q->next; } 

// compariuq chromotJotiuMi and increasing the copies of winner. 

i t ( p :>next ” >revenue>q’' >next - >revenue) p- >next- 

>copi en 4 + ; 

else q- >next ” >copies++ ; 

X 4 ♦ ; 

} 

} 

p-head ; X »d,) ; 
q-diead >uext ; 
p~head; 

while { (ql ^Ntn.L) (q ^ncKpopsize) ) 

{ 

i f (q-'" >copieaa«^0) 

{ //removing chromosomes with 0 copies for next 

generation 

p«>next«q->next ; 
free (q) ; 
q»p«>next; 

} 

else 

( 


‘iiext: ; 


P q; 

*1 p 


p^heiHi ; 

wh i U ' (i ’ 


• ! NPPli) 

//nMk.nr, ropien lor chromosomes 
-■ P “''’Pioa; 

i l !M)) 

f (,)! (x- 2 ; X<, z ;x t + ) 


chrom) ) ) 


} 

pKhi*ad; 
wh i I i*. { |) 


q-(8Cruct chrom *) (malloc (sizeof (struct 

q- >revenue=p- >next ->revenue ; 

q->copies=:0; 

f O r ( y 0 ; y < NIJM^WEEKS ; y + + ) 

for (k«o7k<NlIM_SCR;k+ + ) 

[k] =p->next->b[y] [k] ; } 
q *>naKt*p»>next->next; 

P '"aiext >nextwq; 


p-p aiext ; 

X * 0 ; 

HWKi. 1 

p " - >uu«x; 

p ; 

pi*p ’ >ru!xt ; 

X ♦ ^ ; 


p»ha<id; 

) 


//evaluation function evaluates revenue corresponding to each 
schedu 1 

float evaluationdnt b tNUM WEEKS] [NUM SCR] ) 

( 

int X, y , z , k, xx, yy , run_len9t:hptemp,xxx; 

float revenue^»^0 , 0; 

for { X 1) ; k<NUM weeks ; x 4 1 ) 

for (y- 0 ; y^ NPM SC’H ; y+ + ) 

{ t\in UMiqth«^0; 

k^btx) lyj ; // finding movie which is scheduled in b[x] [yl 
for ( XX«0 / XK<X ; XX-t 4 ) 

foi (yy«0;yy^.y/yy4 O 

it (blxKl tyyl run_length'f+;// calculating run length 

for movia 

run 1 ength 4 ; if (k<6 ) run_length++ ; 

// for firet aix movioB they are one week older than the first week 
so their runltniqth is increased by 1 . 

XXX -X ''rel^dt [k) ; 

//comparing demand of movie with screen capacity. 

i 1 (demand [movie^type [k] ] [xxx] > screen_cap [y] ) 
temp- tjcreen_cap[y] ; 



t ciuis riniJuuMi Ininvie t:ype[k]] [xxx]- 

//ca 1 t'\i 1 ill 1 iMj ifvrana*. 

^ 1 i-vinuH‘ t { ftnp*(’uni ract t:erm[movie 

ret ui u ( i rvrine*) , 


.^ype[k] ] [run_length] ; 


// Imea uiH t iiMt Mt*lec't:y randomly two 

time, peiu'iat ea r^unlum nuuiber; it less than pc, 
site 1 *.intlninl y and t|riHnatey two children, 
void c 1: tn; .nv>ve r ( ) 


chrornosomes each 
selects crossover 


int K, y , xx, z , k, r , r I , r2 , 1 , temp [NUM_WEEKS+io] rMTiM 

struct: chrcm ^ - 

int Cl unu iHAKMI.ind IhOO! . t i tuj, f: 1 agl , f lag2 , m, n • 

int Cl ^ainpvi i i (MAXMl.U'ri 1 SOU] ; ' ' 

peps i zi* i p*lu,vad :uu»xi ; x*0; 
whi lc(pl Nlfbid 


SCR+iO] , count 


X » ♦ ; 

p-p .aicxl ; 

) 

I m ( X 0 ; K* pr:»t»M Izf* i ) cross [x] asx ; 

Id t X ■ U ; X * . I .i| 'iS \Zi* ■ 2 ; X t O 

i 

t - rand ( ) ; 

cnam|<ai t (x! ^>croBs [r] ; 

1 in (y 1 ; y%z ;y+ + ) cross [y] a»cross [y+i] ; 

2 I ; 


an oiispii I i 1 X 1 -ct mis (0) ; 

1*1 oss|fa 4 1 1 X ♦ 1 I -a nmis 1 U ; 
i d' (X- 0 ; X'^ pepsi zv I ; X f- O 

( y^crmispair [x] ; 

p#*lit^4d;q*head/ 

wh i 1 «i ( p - >next - >no 1 «y ) { p»p- >next ;}/ /selected first 

parent 

z f'! rose pa i r [ X 4* 1 1 ; 

wlri ie (c| ' >naxt->nol *2) (q»q->next; }//selected second 


parent 


r rand ( ) I 101;/ /generated random number 


if (r<pc) 

( 

r 1 .Jtid ( ) % MA,X.';i,i)T; // random number for crossover site, 
t (i 1 HiNUM^ iH'H) ; 

1 1 done cruHHovoi between two parents 
1 01 ( k .. t ; k. NUM.. WKKKS ; k + + ) 
tor ( 1 ~0 ; UNtlM SCk;l 4 +) 

temptk) [1] ■•p->next->b[k] [11; 

} 

f Oi ( k t 2 ; k. MUM WEEKS ; k+ + ) 
ioi (l-O; NNUM SCK ;1 + + ) 

p -tiexi >blkl tn-q >next->btk] [IJ ; 


loi (k t WHMK:;; k . t ) 

t Ml t 1 t) 1 • HUH Ml 'K ; 1 . t ) 

( 

i( .iicxi <h i k M I ) 1 Ik! ( 1 ] ; 


XI t ; 

) 

pashetul; 

while (p >t»ext 1 -NULI.) 

( // uueH aame co>. i fi-t i on f unction 

for (xx - (1 ; XX- NIIM KfK 2;xx++) 
for (y»2;y».NUM WKKKy;y* *■) 
for (z-O; Si'sNPM ::C!<,-15 i I ) 

( 


^8 that of initial population. 


{ 1 at 1 1 * 1 ; 

t ’ p ■- >b (y) (*/.] ; 

tin {k-tipk ^ mm ♦) 

I i ii « »t\vni 4i(y 1 1 ik) ) ( f lags«0;flagl»0; 
1 n 1 I i.tq* 1 1 


( 


} 


I or ;iii* y 1 ;mi ♦ ) 

tin {k- (1, k. NUM .SC'H;kt f) 

if (i--p >iie*xt -blml (kl ) { f lagaO;break; } 


1 1 ( (f 1 .UI--U! i.<, ( 1 i.M) 1 - -n 1 
{ 


fount - 0 ; 

itit { 111 t) ; III- NUM WEEKS ; m+ + ) 

tin (k-0;k.N(IM ttf'K.-k+f) 

it (t --p -in-xi .»b[m] [k] ) {count + t; } 

k-O; 

i t (fuum .NUM WEEKS/2) 
whi l u (t IU1.32--O) 


Hd*pl~a; 

n p -ni-xt: >b[y-l) [k] ; 
tin ( 1 .0; UNUM..SCR;1 + +) 

1 1 (n.. -.p >next - >ta [y] [1] ) f lagl=l; 
lid lugU-D) 

p »next. ->b|;y-l] [k] .p->next->b[y] [z] ; 

1 1 a<|2- I ; 


k* * ; 


i‘hie 

whi ie ( i ldq2^^0) 

{ 

t 1 .jg 1 0 ; 

u-p -luixt i*bly l l tk] ; 
tin ( l^U; l^mm 

if {ll- p ,aH»xt *b[yl ll])fi^9i®it 

if ( { 1 at j 1 V ^ 0 ) 



p- >next->b [y] [z] =n; 
f laq2~ 1 ; 

} 

k f + ; 


} 

p.“p- >next / 

} 

p=head; 

while (p- >next 1 =:NULL) 

{ 

p - >next->revenue ^evaluation (p->next->b) ; 
* >next ; 

} 

} 


//mutation function, for each bit, generates random number, checks 
if<pm* if yes, replaces the movie at that bit with movie which was 
not at all scheduled. If no such movie present then selects a movie 
which is released be not played on some other screen in the same 
week. If random no, la <pm/2, screens between movies of same week are 
exchangtKi . 
void mutationO 
{ 

struct chrom *p, *q; 

Int 

x,y, z, k, r,m,n, 1, flag, f lagl, flag2,xx,yy,ch_mov,rr,templ,xxx,count; 
p^head; 

while (p->nextl«NULL) 

{ 

for (xx=0;xx<NUN[_WEEKS;xx+ + ) 

{ 

for (yy=0 ;yy<NUM_SCR;yy++) 

{ 

r=rand() %101;/ /random no. 
if (r<pm) 

for (x=0 ; x<NUM_MOVIES ;x++) 

{ 

f lag«0; 

for (y*0 ;y<NUM_WEEKS ;y++) 
for (Z*0;Z<NUM_SCR;Z++) 

if (x:^=p->next“>b[y] [z] ) 

( 

flag=l; 

} 

} 

if ( (f lag= = 0) && (rel_dt [x] <=xx) ) {ch_mov=x; break; } 
if (flag=xO) p->next->b [xx] [yy] =ch_mov; 


// replaced by movie which was not at all scheduled 

e I « e 

(//replaces by movie which was not scheduled in this week 
fui (x (l;x<NIIM MOV[E3;x(+) 

{ 

1 1 aq 0; 

lot (z ();z<NUM SCR;z + + ) 

( 

if (x===p->next“>b[xx] [z] ) 

( 

flag=l; 

) 

} 

if ( ( f lag=-0 ) && ( rel_dt [x] <=xx) ) (P'->next->b [xx] [yy]=x; break; } 

) 

i f ( r<pm/2 ) 

{ // exchanging screens. 

rr=rand ( ) ‘iNUM^SCR; 
templ« p”“>next->b [xx] [rr] ; 
if ( (templ>NUM MOVIES) | | (tempKO) ) 

{ 

pr intf ( **oh no" ) ; 
getx:h ( ) ; 

) 

P“->next:->b [xx] [rr]“=p->next->b[xx] [yy]; 

if ( (p >next->b[xxl [rr ] >NUM_MOVIES) | | (p->next->b [xx] [rrl<0)) 

{ 

printf(”oh noX"); 
getchO ; 

) 

p->next“>b [xx] [yy]«templ; 

if ( (p- >next '’‘>b [xx] [yy]>NUM MOVIES) | | (p->next-’>b[xx] [yyl<0) ) 

{ 

printf(''oh no2”) ; 
getch ( ) ; 

1 

) 

if ( (gen==®maxgen-l) ) 

( 

print f ( "\n" ) ; 
for(y«0;y<NUM WEEKS;y++) 

f or (2=0; z<NUM SCR;2t+) 

if ( (p">next->b[y] [ z ] >NUM_MOVIES ) 1 | (p->next->b(y] [z]<0) ) 

printf("\ni am hereid no. %d'',p->next->b[y] [z] ,p->next->no) , 
getch { ) ; 

) 

) 

) 

) 

) 

} 

//using same repairing function 
for (xxx=0 ; xxx<NUM_SCR-2 ; xxx++) 



f o r ( y - 2 ; y < NUM _WE E KS ; y (• + ) 
for (z=0; z<NUM_SCR; z I +) 

( 

flag 1= 1 ; 

r = p >nexL-->b[y] [zj; 
flag 1 ; 

for ( k-0; k<NUM_SCR; k h f) 

if (r=-p- >nexL~ >b[y-i] [k] ) ( fldg=0;f lagl=0; } 
if (f lag==l) 

{ 

for (m=0;m<y-l;m++) 

for { k^-0; k<NUM_ SCR; k+ + ) 

if (r==p->next->b[m] [k]) {f lag=0; break; ) 

) 

if ( (flag==0) && (flagl==l) ) 

{ 

count==0; 

f o r ( 0 ; 1 u < y - 1 ; m + + ) 
for ( k^O; k<NUM SCR; k+f ) 
it >next->b[in) [k] ) {count++; ) 

k-0; 

flag2==^0; 

if (count>NUM_WEEKS/2) 
while ( f lag2===0) 

{ 

flagl^O; 

n=p“>next->b [y-1] [k]; 
for (1-0;1<NUM__SCR;1++) 
if (n==p->next*“>b [ y] [1] ) f lagl=l; 
if (flagl«=0) 

I 

P“>next->b [y-1] [k] ~p~>next->b[y] [z]; 
flag2=l; 

) 

k++ ; 

) 

elae 

while (f lag2«“0) 

I 

flagl^O; 

n«p->next“>b[y-l] [k]; 
for (1«0;1<NUM_SCR;1++) 

( n=«*p-'>next->b [y] [1] ) f lagl=l; 
if (f lagl*“0) 

{ 

p->next“'>b [y] [z]~n; 
f lag2«l; 

! 

k-H^; 

) 

1 

1 

if ( (gen*«maxgen-l) ) 

{printf ( ”\n” ) ; 

for (y=0;y<NUM_WEEKS;y+f ) 

{ for (z«0;z<NUM_SCR;z++) 



ii( (p-next >b[yl [ 2 ] >NUM_MOVIES ) 1 1 (p->next- >b [y ] [2]<0)) 

printf ("\n'i,d no. ‘id" , p->next->b [y] [ 2 ] , p->next->no) ; 
getch ( ) ; 

) 

) 

) 

p=p->next ; 

) 

p=head; 

while (p->next !=NULL) 

( 

p->next->revenue=evaluation(p->next->b) ;p=p->next; 

) 

) 


Sample data file giving the type and release date for movies: 

21113 1 . 11433134213414313312414313132 
3 13 

- 1 - 1-- 1 - 1 - 1-1000011112 22233334 4 4455556 
6 6 6 7 7 7 7 
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